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1.  Summary. 


This  report  considers  the  use  of  higher-order  statistics  (HOS),  in  particular  the  triple-correlation 
function  (TCF)  and  ninth-order  correlation  function  (NCF),  to  detect  the  presence  of  m-sequences 
and  Gold  codes  in  signal  intercepts.  Section  2  describes  the  continuation  of  work  in  [1]  to  blindly 
detect  m-sequences  from  the  TCF.  The  new  work  reported  in  2.4  improves  detection  rates 
significantly  by  searching  for  the  TCF  patterns  associated  with  m-sequences  likely  to  be  present. 
Section  3  presents  a  new  theory  for  the  TCF  of  Gold  codes.  The  expected  TCF  of  Gold  codes  is  3- 
valued:  no  peaks  of  constituent  m-sequences  are  present.  The  detection  of  a  constituent  m-sequence 
of  a  Gold  code  relies  on  the  property  that,  at  the  locations  of  that  m-sequence’ s  TCF  peaks,  the 
Gold  code  TCF  always  has  the  expected  value  of  -1/L,  L  the  sequence  length.  Thus  the  Gold  code 
TCF  may  be  searched  for  sets  of  these  values  at  known  TCF  peak  locations  for  likely  m-sequences. 
When  a  constituent  m-sequence  is  found,  a  preferred-complementary  m-sequence  and  the  actual 
Gold  code  may  also  be  easily  found.  Section  4  examines  higher-order  (>3)  correlation  functions.  It 
is  shown  that  Gold  codes  have  NCF  peaks.  Although  the  NCF  is  very  noisy,  there  are  many  peak 
locations,  derivable  from  the  TCF  peak  locations  of  the  preferred  pair  of  m-sequences,  over  which 
to  average.  The  NCF  method  has  the  potential  to  produce  good  blind  detection  at  high  noise  levels 
for  long  Gold  codes.  The  final  methods  of  section  5  were  suggested  by  R  Gold  after  a  presentation 
by  the  author  at  AFRL,  WPAFB  in  December  2002:  a  1 -stage  triple-correlator  detector  for 
m-sequences  and  a  2-stage  triple-correlator  detector  for  Gold  codes.  Both  detectors  may  be  simply 
implemented  with  delay  circuits,  but  their  outputs  are  extremely  noisy.  They  also  provide  the  basis 
of  phase-independent  decoders. 

With  the  exception  of  Dr  Gold’s  suggestion  and  the  simulation  results  for  m-sequences  in  2.3, 
produced  by  K  Batty,  all  material  reported  is  the  original  work  of  the  author. 


2.  Detection  and  Identification  of  m-sequences. 

It  has  been  shown  that  m-sequences  in  signal  intercepts  may  be  detected  and  identified  using 
HOS  [1].  Subsequent  work  has  demonstrated  it  is  possible  to  automate  the  blind  (unknown  length 
and  possible  m-sequences  present)  derivation  of  an  m-sequence’ s  generator  polynomial  from  a 
single  TCF  peak  (see  2.1).  Extensive  simulation  has  revealed  how  intercept  length  and  noise  level 
influence  detection  rate  (2.3).  This  extensive  simulation  was  made  possible  by  assuming  hard 
decoding  by  the  receiver:  the  chip  values  processed  were  +1  or  -1  rather  than  real  numbers, 
allowing  fast  binary  arithmetic  (2.2). 

The  simulation  results  show  that  the  blind  algorithms  are  inadequate  for  LPI  signals,  i.e. 
identification  rates  are  low  for  SNRs  below  -5  dB,  unless  long  or  multiple  intercepts  are  available. 
The  present  research  is  investigating  the  alternative  strategy  of  searching  for  known  m-sequences 
with  known  TCF  peak  locations  (2.4). 


2. 1  Automatic  blind  m-sequence  identification. 

Firstly,  the  partial  TCF  C'(p,q )  of  an  (V-length  intercept  x,  is  calculated, 

j  N~‘i 

c'(p>  <?)  =  E  xixi+Pxi+q  ’  0<p<^<A-l, 


C'(^,p)  =  C'(/?,^)  for  0<q<  p  <  N  -l. 


(2.1) 


and  peaks  are  detected.  Each  peak  is  used  to  estimate  the  length  of  an  m-sequence  present, 
and  a  majority  vote  taken  to  determine  the  final  sequence  length  L.  A  single  peak  is  then 
selected  to  derive  the  generator  polynomial  g(X)  of  the  m-sequence;  other  peaks  may  then 
be  used  to  verify  g(X).  A  peak  location  (p,q )  allows  the  algorithm  to  determine  a  unique 
g(X)  if  the  following  criterion  is  satisfied: 

gcd(p,L)  =  1  v  gcd(?,L)  =  l,  (2.2) 

where  gcd  is  greatest  common  divisor  and  v  is  logical  ‘or’.  Note,  when  L  is  prime,  all  (p,q ) 
determine  a  g(A). 

Two  algorithms  have  been  developed  [2]  to  derive  g(X)  from  a  (p,q )  satisfying  (2.2). 

(i)  Gauss  elimination  over  GF(2)  for  the  m  equations: 


er  +  aq  =  1, 

^,2/>inodL  £^2(/mod  L  _ j 


a 


2m~'  p  modi  q2"'  't/mod  L  _  j 


(ii)  Use  the  derived  equation: 


(2.3) 


ap  +  aq  - 1  =>  ap~q  +  a~q  =  1  (2.4) 

to  rewrite  any  ot,  in  equations  derived  from  peaks,  in  terms  of  reduced  powers  of  oc. 

a"  =an~q(ap  +l).  (2.5) 

Fault  tolerant  software,  using  either  (i)  or  (ii),  has  been  written  which  has  produced  a  g(X)  in 
all  cases  tried. 

More  details  of  the  above  methods  may  be  found  in  [2,3]. 


2.2  Noise  assumptions  for  hard-decision  receivers. 

HOS  suppress  the  influence  of  Gaussian  noise  signals,  their  cumulants  above  order  two 
vanishing.  A  partial  TCF  is  a  third-order  moment  which  closely  approximates  the  third- 
order  cumulant  as  the  expected  average  value  of  a  partial  m-sequence  is  -1/L,  close  to  zero 
for  large  L.  For  processes  with  zero  mean,  third-order  moments  and  cumulants  are  identical. 
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The  statistical  properties  of  partial  TCF  estimators  are  derived  in  [1]  for  m-sequences  whose 
values  are  corrupted  by  AWGN.  For  N(0,CT)  noise,  the  additional  variance  of  C'(p,q )  is 


var  , 

O' 


V+3<74+3cr2N 

v  1  , 


for 


p±q. 


(2.6) 


An  alternative  noise  model  results  if  the  receiver  makes  a  hard  decision  for  each  chip, 
assuming  the  same  AWGN.  Sufficient  noise  may  result  in  a  transmitted  1  being  received  as 
-1  or  vice  versa.  Assuming  a  threshold  of  zero  for  the  hard  decision,  the  probability  of 
mistaking  1  for  -1 ,  or  -1  for  1 ,  is 


P 


dr. 


(2.7) 


Thus  the  3-valued  chip  error  y,-  has  the  probability  distribution: 

dy,-=o)=i-2# 

P(y,=  2)=fi,  (2.8) 

P(y,=-  2)  =  A 

The  partial  TCF  of  x,  =  w,  +  y, ,  for  m-sequence  sample  u,  with  superimposed  error  sequence 

y; ,  is 

c'.v (p,q) = -rr~  i! («,•  +  >v  tui+P  +  y>+P  Xui+q  +  y,+q ) 

N-q  i= i 

=  C’„(p,<l)+C'y(p,q)+iTJ  (2.9) 

j= 1 


Cases  p  =  q,p  =  0  and  q  =  0  are  not  considered,  as  there  are  no  peaks  for  such  values. 


C'u(p,q):  the  partial  TCF  of  u 

is  unbiased,  with  variance 

var  [C'tl{p,q)]  =  < 

1 

N-q 

for  non  peak  ( p,q ) 

(2.10) 

0 

for  peak  ( p,q ) 

C'y{p,q ): 

eLv,.>',+;,v)+J  = 

0, 

var[yj  =  8/?. 

(2.11) 

Thus 

e[c;  (/>,«)]=«, 

var [c  {p,q)\=  &  . 

}  N-q 

(2.12) 
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Ty  T2,  Ty 

T\i  =  Uj+k  yi+q 

for 

Mj+k  —  MjUj+p  , 

Ty  —  Mj+iVj+p 

for 

Uj+l  —  UjUj+q  , 

Ty  =  Uj+Syi 

for 

Mj+s  —  Llj+p  Uj+q. 

T],  T2  and  T$  are  nearly  orthogonal  combinations  of  y  values: 
E[r„]  =  0 


var[r„]  =  E[u’11  v;li;]=R/? 

Clearly  Ti,  and  Ty  have  the  same  mean  and  variance,  and 


eM=o 


yi  <  i  <  3 


var[7)]  = 


_8^_ 

N-q 


(2.13) 


(2.14) 


(2.15) 


T4,  T? Te'. 

T\  i  —  Uj  yj+p  yi+q  5 

T$j  —  Uj+p  V/  yi+q  5 

%i  =  ul+qytyi+p  .  (2.16) 

Again,  the  7}  are  nearly  orthogonal  combinations  of  the  products  of  independent  y  values: 
Efcj  =  0, 

(2.17) 

var[r4,  ]  =  E [uf  yf+p  yf+q  ]  =  64/?2  , 


so 

E[rJ=E[r5]=E[7-6]=o, 


vaifc  ]  =  var[l5  ]  =  varfc  ]  =  . 

N-q 


(2.18) 
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C[  v  is  therefore  an  unbiased  estimator  for  w’s  TCF  for  pt-q,  p^O  and  q£ 0.  As  the  above 

contributions  are  nearly  independent,  the  overall  variance  for  non-peak  (p,q)  is 
approximately: 


(2.19) 


For  peak  (p,g)  this  variance  (2.19)  is  reduced  by  1  /(N-q).  These  approximate  variances 
were  used  in  algorithms  to  detect  peaks  and  determine  sequence  length  L. 

TCF  computations  are  significantly  faster  when  noise  is  represented  as  chip  errors. 
Multiplications  may  be  replaced  with  XORs  of  the  binary  data  in  the  C  programs. 


2.3  Simulation  results  for  blind  identification. 

A  set  of  simulations  was  conducted  to  investigate  the  joint  effects  of  probability  of  chip 
error  P(e)  and  intercept  length  N  on  identification  rate.  The  m-sequences  used  were  of 
length  L-  127,  knowledge  denied  to  the  algorithm.  A  series  of  1000  trials  was  carried  out 
for  a  set  of  P(e)  varying  from  0.01  to  0.2  and  N  varying  from  60  to  1 20.  Results  are  shown 
in  Figure  2.1. 


Number  of 
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1000 

900 


100 


p(Chip  Error) 


100 
B0 

60 
OJ 

o 


■  900- 1000 ! 

■  800-900  j 
0700  800  j 
|  ■  600-700  | 
|  Q  500-600  [ 

■  400-500  [ 
In 300-400  i 
!  0200-300  I 

■  100  200  : 
no  too  j 


120  Intercept 
Length 


Figure  2.1  Influence  of  chip  errors  and  intercept  length  on  blind  g(X)  identification. 


Note  that  P(e)  =  0.06  corresponds  with  a  SNR  of  approximately  3  dB.  For  3  dB:  when 
N  =  1 20  the  identification  rate  is  96%;  when  intercept  length  is  halved  to  N  =  60,  the 
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identification  rate  is  almost  halved  to  5 1  %.  Even  for  large  N,  close  to  L,  identification  rates 
are  poor  for  large  P(e),  or  low  SNR.  They  may  only  be  improved  by  using  multiple 
intercepts. 


2.4  Pattern  recognition  to  identify  known  m-sequences 

The  number  of  primitive  polynomials  S„,  of  order  m,  each  generating  an  m-sequence, 
satisfies: 

2'"  -  2 

Sm  = -  when  L  =  2"'  - 1  is  prime, 

m 

(2.20) 

Sm  <  —  for  non  -  prime  L. 
m 

Very  few  of  the  TCF  peaks  of  the  generated  m-sequences  are  common  with  those  of 
m-sequences  of  the  same  or  different  length. 


The  total  number  of  primitive  polynomials  up  to  order  13,  i.e.  m-sequences  up  to  length 
8,191,  is  1,108: 

5X  =U08.  (2.21) 

m 

This  is  not  a  large  search  space.  The  number  of  primitive  polynomials  up  to  order  of  17  is 
7,710;  the  number  up  to  order  19  is  27,594. 


Discriminant  functions  of  the  form 


D.=tlX,C’(p„ql)  (2.22) 

(=1 

were  derived  for  a  set  of  m-sequences  of  length  1 ,023.  Each  of  the  nu  peaks  (p,-,g,- )  are 
unique  to  sequence  u.  Peaks  close  to  the  origin  were  chosen  so  they  could  be  evaluated 
from  short  samples.  The  weightings  A  (2.22)  of  TCF  values  at  peak  locations  i  reflect  the 
variance  of  those  values  {approximately  l/(/V—  q,)}  and  modulation  effects. 


Simulations  to  estimate  identification  rate  were  carried  out  for  noisy  m-sequences  and  noisy 
modulated  m-sequences.  For  unmodulated  sequences,  the  A  were  inversely  proportional  to 
variances: 

A i=N- qt  ,  q,  >  Pi .  (2.23) 

For  randomly  modulated  sequences,  the  A  were  estimated  by  multivariate  analysis.  In  each 
run,  the  sequence  was  identified  when  the  appropriate  Du  (2.22)  exceeded  a  threshold.  Two 
intercept  lengths  were  used:  127,  when  the  Du  were  linear  combinations  of  5  peaks  ( nu  =  5); 
255,  when  20  peaks  were  used.  The  results,  summarized  in  tables  2.1  and  2.2,  are 
significantly  better  than  those  for  blind  identification  in  section  2.3. 
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SAMPLE  / 1023 

127  (5  peaks) 

255  (20  peaks) 

SNR 

3dB 

96.7  (0.5) 

100(0) 

OdB 

88.3  (1.3) 

99.3  (0) 

Table  2.1  Identification  rates  for  m-sequences  with  no  modulation. 


SAMPLE  / 1023 

1 27  (5  peaks) 

255  (20  peaks) 

SNR 

3dB 

84.5  (0.7) 

99.5  (0) 

OdB 

76.1  (1.2) 

97.4  (0) 

Table  2.2  Identification  rates  for  randomly  modulated  m-sequences. 

For  the  tabulated  x(y)  values,  x  is  the  %  identification  of  an  m-sequence  and  y  the  %  false 
alarms,  the  detection  of  a  sequence  not  present. 


3.  Triple-Correlation  Based  Detection  and  Identification  of  Gold  Codes 

Gold  codes  are  more  difficult  to  detect  than  m-sequences  from  their  triple-correlation  function. 

They  are  not  finite-field  elements  and  do  not  obey  the  shift-and-multiply  rule.  There  are  no  triple¬ 
correlation  peaks  except  when  sequence  length  L  is  divisible  by  3,  when  the  two  peaks  (L/3,  2L/3) 
and  (2L/3,  L/3)  on  the  diagonal  p  +  q  =  L  are  shared  by  all  L-length  sequences,  e.g.  for  L  =  15,  there 
are  peaks  at  (5,10)  and  (10,5). 

Gold  codes  have  a  looser  structure  than  m-sequences:  shift-and-multiplying  produces  not  another 
cyclically  shifted  version  of  the  same  code  but  another  code  in  the  same  family,  i.e.  one  generated 
by  the  same  two  preferred  m-sequences  but  with  a  different  relative  phase  (see  3.2).  It  will  be 
shown  that  the  triple-correlation  function  of  Gold  codes  contains  sufficient  information  for  the 
detection  of  constituent  m-sequences  and,  by  further  correlation  processing,  the  actual  Gold  code. 


3.1.  Preferred  m-sequences. 

Gold  codes  are  formed  by  modulo-2  adding  the  bits  of  a  preferred  pair  of  m-sequences. 
Preferred  pairs  have  the  defining  property  that  their  cross-correlation  function  is  3-valued. 
For  a  preferred  pair  of  m-sequences  of  length  L-2"  - 1 ,  these  three  values  are  [4]: 

0_x  =  —t(n)/L 

0O=~\/L 

9\  ={t(n)-2}/  L  (3.1) 
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where 


t  («)  = ' 


1  +  2 
1  +  2 


(ii+l)/2 

(n+2)l2 


for  odd  n 
for  even  n 


(3.2) 


Note  that  there  are  no  preferred  pairs  for  n  =  0  mod  4. 

Preferred  pairs  may  be  constructed  by  decimation.  For  sequences  u  =  (u0,  uv  ..,  uL_ ,)  and 
v  =  (v0,  v„  vw). 


v  =  u[q ]  =>  V/(0  <  i  <  L - 1)  v,  =u 


qi  mod  L  * 


(3.3) 


u[q ]  is  a  proper  decimation  of  u,  i.e.  has  period  L,  iff  gcd  (L,  q)  =  1 . 
w  and  v  are  a  preferred  pair  if  n  ^  0  mod  4  and  v  =  w[g],  where  q  is  odd  and 
4  =  2*+!  or  q  =  22k  -2*+l„ 


gcd  (n,  /c)  = 


1  for  odd  n 

2  for  n  =  2  mod  4 


(3.4) 


For  L  =  3 1  (n  =  5),  the  preferred  pairs  of  m-sequences  are  those  connected  by  a  straight  line 
in  Figure  3.1  [4]. 


w[15] 


Figure  3.1  Preferred  pairs  of  m-sequences  for  L  =  3 1 . 


For  u,  45  is  the  octal  form  of  its  generator  polynomial’s  coefficients: 


g(,(X)  =  X5+X2+l 


(3.5) 


3.2.  Derivation  of  Gold  codes  from  preferred  pairs  of  m-sequences. 

Gold  codes  may  be  generated  from  a  pair  of  preferred  sequences  u  and  v  by  modulo-2 
adding  u  and  cyclically  shifted  versions  of  v,  or  vice  versa.  If  T  denotes  the  operator  which 
cyclically  shifts  vectors  one  place  to  the  left: 

7v  =  (vpv2,  ..,  vL_,,v0), 


r*v  =  (v*,v*+1,..,  vL_,,v0,..,  vt_,)  for  0  <k<L 


(3.6) 


then  the  following  set  of  L+2  Gold  codes  results: 

G(w,v)  =  {u,v,  u  ©  v,  u  ©  Tv,  u  ®T2v,  ..,  u  0  TL~]v  }  (3.7) 

where  0  denotes  modulo-2  addition  of  corresponding  sequence  bits.  For  the  L  Gold 
sequences 

wk  =  u  ®Tkv,  0<  k  <  L  —  \  (3.8) 

individual  bits  may  be  expressed: 

=  M,  ©  ^(/+A)mod L’  0  <  /  <  L  -  1  (3.9) 

As  an  example,  the  preferred  pair  of  sequences  u  and  v  =  w[3]  in  Figure  3.1 ,  octal  45  and  75, 
and  vv°  =  w©v  are: 


u= 1 0 1 0 1 1 101 10001 1 1 1 1001 1010010000 
V  =  n[3]  =  1011010100011101111100100110000 

w°  =  0001  101  1  1  101  1010001  1  1  1  1  10100000  (3.10) 


Note  that,  if  bits  are  mapped  {0 1 — >  1, 1 1— >  —  1 },  the  ±1  values  are  multiplied: 


W:  =U;Vt 


(i+k) mod  L  • 


0  <  /  <  L  - 1 


(3.11) 


e.g. 

w°  =  11 1-1-1 1-1-1-1-11-1-1 1-1 1 1 1—1—1— 1—1— 1—1 1-111111 


(3.12) 
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3.3.  Autocorrelation  functions  (ACFs)  of  Gold  sequences. 

For  the  {1,-1}  representation,  assuming  all  indices  are  mod  L,  the  ACF  of  wk  is  given  by 


c‘(P)4i>.v 

L  I- 1 


i+P 


i=0 


]_ 

L 


L-\ 

I«,h+t«i+/,+hp 


i=0 


(3.13) 


1  L~ 1 

7  ZL  ^'  +  'V  ^  ,+  Vl/ 
L  /=0 


for  some  r  e  {0,l,..,L-l},  rp  determined  by p  and  ,va7,  by  k  and p.  Thus  (4(p)  reduces  to 
the  cross-correlation  between  two  preferred  m-sequences,  and  has  the  values: 

(4(0)= 1 


Ck  (p)  =  0_x , 6q  or 6] ,  as  given  in  (3.1),  for  1  <  p<  L- 1 


(3.14) 


Partial  correlation  functions  of  Gold  sequences  may  take  a  wider  range  of  values,  which 
depend  on  initial  position  in  the  sequence.  Thus  the  following  partial  ACF  of  wk,  calculated 
from  an  A-sum  ( N<L ),  depends  on  j  and  p: 


Li  i=j 


<+/' 


(3.15) 


But,  using  group  closure  (shift-and-multiply  rule)  for  m-sequences  u  and  v: 


wf  wlP  =UiVl+kUl+pVM+p 


=uiui+pvi+kvi+k+p 

=ui+rvi+hv,  rp  and  skp  e  {0,1,  ..,L-l} 


(3.16) 


Thus 


CkN{j,  0)  =  1 

(3.17) 

i  N+j- 1 

cUj.P)=-’Zwk- 

iV  i=j 

The  partial  correlation  of  wk  reduces  to  the  mean  of  a  partial  sequence  within  another 
v4from  G(u,v). 
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3.4.  Triple-Correlation  Functions  (TCFs)  of  Gold  Sequences. 


3.4.1  TCFs  of  Complete  Gold  Sequences 

The  TCF  of  wk  is  defined: 


1  L~' 

Ck{p,q)  =  -Yuw"wlvwlq 

L  i= o 


(3.18) 


Substituting  from  (3.16) 


1 


Ck(p,q)  = 

L  ,= o 


i+q 


(3.19) 


The  TCF  is  thus  the  cross-correlation  of  two  Gold  sequences  from  G(u,v).  This 
cross-correlation  is  3-valued  [4],  with  the  values  given  in  (3.1),  so 

Ck(p,q)  =  0_],0oor0]  (3.20) 

(3.1 8)  may  be  rewritten: 

Ck  (p,q)  =  jY,uivi+kU'+rvM+nui+‘,vi+k+<i  (3-21) 

L  i=0 

Consider  the  case  where  ( p',q' )  is  a  TCF  peak  location  for  m-sequence  u, 
i.e.  Uj  =  ul+p.ut+q.  for  0  <  i  <  L  - 1 .  As  u,ui+p. ui+q.  =  uj  =  1 : 


1  i_1 

Ck{p',q)  = 

^  i=0 


(3.22) 


For  prime  L  and  nearly  all  other  cases  when  L  is  not  prime  [1],  if  ( p',q ')  is  a  TCF 
peak  for  u  it  cannot  be  a  TCF  peak  for  v,  i.e.  it  is  impossible  that 

viu  =  vM+p-vMW  for  all  i  (0  <  /  <  L  - 1)  (3-23) 

Closure  (shift-and-multiply  rule)  gives 

Vi+k+P'  Vi+Uq  =  Vi+k+r,lV  .  V?  ^  0  (3’24) 

for  some  r  - .  (1  <  r  •  .  <L- 1)  determined  by  (p',q') . 


ll 


Thus,  as  phase  is  irrelevant  to  the  TCF: 
Ck{p',q')=jYJvi+kvi+k+ri., 

L  i'=0 


1  L~] 

= — S"1  vv 

J  La  I  '+V 

L  i=0 


(3.25) 


But  this  expression  for  the  ACF  of  v  yields  just  one  value  as  rp.q.  cannot  be  zero: 


V,-*0=»C*(/>V)=y  (3-26) 

A  similar  reduction  occurs  when  ( p',q ")  is  a  TCF  peak  for  v,  and  automatically  not 
for  u: 

Ck{p\q")  =  \YJuiui+pMi+q. 

L  i= 0 


1  L~l 

=~y ,  s..  ^ o 

7  l  l+spy  7  p  q 

i=0 


(3.27) 


-1 

L 


Note  that,  from  (3.21)  when  p=q: 


Ck{p,p ) 


1  1-1 
L  1=0 


1  z'_l 

t2>‘ 

^  /=o 


(3.28) 


for  all  p  (0  <  p  <  L- 1) .  This,  the  cross-correlation  of  two  preferred  m-sequences,  is 
3-valued:  6L,,  #0  or  ft  are  the  possible  values  for  different  k. 


Also,  as,  for  all  a,  v2a-u2a=  1, 


1  L~] 

Ck{0,q)  =  -YJui+<,vi+k+<l 

L  1=0 

<3-29> 

^  1=0 

C(0.0)  =  lg«iv,,t 
L  ,=o 
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But  these  three  values  are  equal: 

Ck  (0  ,q)=Ck  (p,0)=Ck  (0,0  )  =  CM{k)  (3.30) 

where  CUl  v)(k)  is  the  cross-correlation  between  u  and  v  for  shift  k,  known  to  be 
3 -valued. 

Thus,  when  p  =  q  or  either  or  both  of  p  and  q  are  zero,  the  same  value  $.u  or 
results,  depending  on  k\ 

p  =  q  v  p  =  0  v  q  =  0  =>  Ck {p,q)  =  C(llY){k)  =  (3.31) 

L  i=0 

The  L  - 1  (p\ q)  and  L  - 1  {p",q")  values  resulting  in  TCF  peaks  for  exclusively  u 
or  v  produce  2L-2  values  of  60  for  the  TCF  of  wk: 

{(p.9)=(p,'.9,')®  (p>9)  =  (p,W)}>  1  <i<L-\  =>  Ck{p,q)=  dQ  =  y- 

(3.32) 

where  ©  denotes  ‘exclusive  or’.  For  all  other  Ck  ( p,q ),  which  may  be  regarded  as  a 
cross-correlation  function  of  two  Gold  sequences  from  G(w,v): 

{p,  q)  *  {p'n  q'i ) A  (p.  <i )  *  (pf»  9,') A  p  *  q  =>  (p.  9)  =  £-1 .  or  3  • 


(3.33) 


If  there  are  equal  numbers  of  (p,q)  values  producing  0_,  and  ,  the  mean  of  the 
resulting  C*( p,q )  is 

(<9_,  +0,)  _ -r(«)+r(«)-2  _ -1 

2  2 L  ~  L 
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Equivalent  forms 


7 1  "L-wL, 

Lt  i 


0-\,  6b  or  Q\ 


Y  ui  vi+fc  Uj+ p  Vj+k + p  Uj+q  V,+ j.  +Cj 
L>  i 


(p>q)=(pW) 

a  TCF  peak  location  for  u , 


i.e.  u-,ui+p'Ui+q'=\. 


j  Y  Vj+kVj+k+p'Vj+k+q' 

L,  i 


(p\  q  )  a  peak  for  u 
=>  NOT  a  peak  for  v 


TY.VjVj+r(r±0) 


L  j 


Oq  - 


-1 


Figure  3.1  Alternative  TCF  Values 


3.4.2  TCFs  of  Partial  Gold  Sequences 

The  partial  TCF  of  wk  may  be  defined  (N  <  L): 


]  N +./-1 

CkN(j,p,q)=—  'ZW'W'+i>W'+<,  ’ 

N  i=j 


which  may  again  be  simplified  using  (3.16): 

1  N+j-] 

CkN  {j,  p,q)  =  —  y  WLP  WL, 

t*  i=j 


(3.35) 


(3.36) 


which  is  the  partial  cross-correlation  between  two  subsequences  of  different  Gold 
sequences  from  G(u,v).  (3.36)  may  be  rewritten: 

|  N+j- 1 

ClN(j,p,q)  =  —  'YuUivMui+pvi+l+pui+qvi+k+q  (3.37) 

tv  i=j 


Rewriting  (3.37)  for  the  general  case,  where  (p,q)  is  not  necessarily  a  peak  for  u  or  v: 

j  N+j- 1 

ckN  ( j ,  p>q)=—  Z  uiui+pui+9vMvi+k+pvM+q 

tV  i—j 


N+j- 1 


=-  y 


N 


Ui+am  Vi+Pkr 


i=J 


(3.38) 


j  N+i i-1 

=  T7  Z  w>  >*  =  j  +  apq  and  m  =  pkpq  -  apq 
tv  ,=/> 

Thus  the  general  TCF  may  be  regarded  as  a  cross-correlation  between  two  Gold 
sequences,  a  cross-correlation  of  a  pair  of  preferred  m-sequences,  or  the  average 
value  of  a  Gold  sequence.  For  a  complete  sequence  (N  =  L),  all  these  are  3-valued. 

For  ( p',q ')  corresponding  to  a  TCF  peak  for  u  and  not  for  v: 

1  N+j-i 

C  kN  ( j ,  p',  q)  =  —  y  v,+,  vi+k+p’Vi+kW 

tV  i=j 


]_ 

N 


N+ j+k- 1 

SVV-  -V-  - 
/  i+/)  i+q 


i=j+k 


j  /V+y+Jt-l 

=  T7  2>.v,.w  ■  r„*0 

i-j+k 


(3.39) 
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Thus  C'N(j,  p',q)  is  a  partial  ACF  value  (for  non-zero  shift)  of  an  m-sequence,  with 
a  mean  of  - 1 /L,  and  its  statistical  properties  independent  of  k  and  j.  Similar  values 
arise  if  ( p,q )  is  a  peak  for  v  and  not  for  u. 

Again,  if  p=q  or  either  or  both  of  p  and  q  are  zero,  the  average  value  of  the  partial 
TCF  is  0_,,  60  or  0,,  depending  on  k. 


3.5  Statistical  theory  of  Gold  sequence  TCFs. 

For  (jj,q )  exclusively  a  peak  for  u  or  v,  the  partial  TCF  in  (3.39)  reduces  to  a  partial  ACF  of 
an  m-sequence.  [1]  gives  the  following  mean  and  variance: 


E [C*  (;>,<?)]  =  y 

(3.40) 

var[c*  {j,  p,  q)]  =  ~ |l  -  ~  j} 

As  an  AMength  intercept  is  available,  N'  =  N-q  products  are  used  to  estimate  CkN  ( j ,  p,  q) , 
assuming  p  <  q. 


For  the  more  general  partial  TCF  expression  (3.38)  when  (p,q )  is  not  necessarily  a  u  or  v 
peak,  there  are  three  possible  cases:  the  mean  value  (actual  value  when  N=L )  may  be  0_,,  0O 
or  d\.  Assuming  each  chip  of  the  Gold  sequence  may  be  regarded  as  independent  of  the 
other  chips,  the  first  two  moments  about  zero  of  the  partial  TCF  are: 


e[cJ  (j,  pA= 


N+A-1 


I 

i=A 


W; 


=  —0(N  =  0. 
N  c 


(3.41) 


fr  N+A-lN+A-l  r 

W?lN(j,p,q)}  \=  2  S  Ekxl 

i=A  j=A 

(3.42) 

= -^r  W + (w2  -  wk  }= ^-{i + (n -!)<?, } 


where 


s,  =E [«,>;] ,  i*j. 


(3.43) 
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So  the  variance  is 


var[c*  (./,  p, g)\  =^{\  +  {N- \)S(  -  0?N }  (3.44) 

Because  this  variance  is  zero  when  N  =  L  (the  partial  TCF  becomes  0, ),  setting  (3.44)  to 
zero  gives 

S  =  Mzi  (3.45) 

f  L- 1 


Substituting  (3.45)  into  (3.44)  gives 


(L-N)ey 

L- i 


N{L-\) 


(3.46) 


assuming  L  is  large  and  6 f  is  small.  The  variance  approximation  in  (3.46)  is  valid  in  all 
three  above  cases.  The  Central  Limit  Theorem  may  be  used  to  show  all  three  populations 
are  approximately  Gaussian  with  the  same  variance  but  means  of  6L],  and  6\. 

For  the  case  t  =  0  {0O=- 1  /L)  the  exact  variance  is: 


L-N  I 


The  expression  in  (3.47)  is  asymptotically  identical  to  the  variance  in  (3.40), 
interpreting  N'  and  N  as  the  number  of  products  used  to  estimate  the  partial  TCF. 


(3.47) 
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3.6  Statistical  tests  for  constituent  m-sequences  of  a  Gold  code. 

In  the  ideal  case  of  a  noiseless  L-length  intercept  of  a  Gold  code,  the  presence  of  constituent 
m-sequences  u  and  v  would  be  indicated  by  the  -1/L  TCF  values  at  all  peak  locations 
( p',q ')  for  u  and  {p',q”)  for  v.  For  an  /V-length  intercept  ( N<L ),  the  mean  TCF  values  at 

peak  locations  is  still  -1/L,  but  with  the  approximate  variance  <702  given  in  (3.46). 

However,  for  a  set  of  TCF  values  from  non-peak  locations,  the  mean  value  would  also  be 
-1/L,  as  the  expected  numbers  with  average  values  of  6-\  and  0\  are  equal,  and  their  average 
is  -1/L  (3.34).  The  approximate  variance  for  such  a  set  is  given  by 


'±i 


=  <Jq+7\ 


t(n)-l 


,L(n+2)/2j  ' 


=  05+ri 


+d2 


(3.48) 


where  LxJ  denotes  the  largest  integer  not  exceeding  x  and  /is  the  proportion  of  TCF  values 
taking  the  values  -t(n)/L  or  {t(n)-2}/L.  The  theoretical  value  for  /is  Vi.  This 
approximation  (3.48)  is  an  undere stimate,  but  is  used  for  convenience  in  the  tests  below, 
which  therefore  yield  pessimistic  detection  rates.  The  actual  variance,  assuming  is  exact, 
is 


Jx2e”('v~d)3/2<T" 


dx 


(3.49) 


It  is  proposed  that  tests  for  u  and  v  be  based  on  the  different  variances  of  sets  of  TCF  values 
at  peak  and  at  non-peak  locations.  For  any  L,  it  is  easy  to  tabulate  peak  locations  for  all  re¬ 
sequences.  The  majority  of  locations  are  non-peak  and  a  set  of  these  may  be  chosen  as  a 
control  set.  An  approximate  test  may  be  carried  out  on  whether  the  variance  of  the  control 
set  is  significantly  greater  than  the  variance  of  the  proposed  set  of  peaks  by  calculating  the 
ratio  of  the  variances  and  comparing  with  an  F-statistic.  The  following  sum  of  squared 
deviations  from  -1/L  is  computed: 


Z  = 


(3.50) 


For  Zu  the  value  of  Z  when  (p, •,</■)  are  /,  peaks  for  u,  and  Zc  when  the  /.  peaks  (pi,qi)  are  the 
non-peaks  in  the  control  set,  the  variance-ratio  may  be  compared  with  tabulated  F  values: 
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0 


TCF  value 


-1 

L 

Figure  3.2  PDF  of  set  of  TCF  values  including  just  peak  locations. 


Figure  3.3  PDF  of  control  set  of  TCF  values  including  no  peaks. 
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(3.51) 


2r/(y,-l) 

z./O'.-i) 


>F, 


(2,-1  )•(/„-!) 


u  present. 


The  test  is  only  approximate  as  the  mixture  of  three  Gaussian  populations  with  different 
means  in  the  control  set  is  not  Gaussian  (see  Figures  3.2  and  3.3  above).  However,  the 
presence  of  AWGN  in  the  original  intercept  results  in  a  closer  approximation  to  Gaussian. 
Note  that  the  test  is  still  valid  if  the  two  populations  of  means  and  6\  are  unequal  as  the 
average  distance  of  any  TCF  value  from  -ML  is  the  same  d  =  { t(/z)  -  1  }/L  in  both  cases. 

As  the  partial  TCF  of  a  Gold  code  at  peaks  for  one  of  the  constituent  m-sequences  reduces 
to  the  partial  TCF  of  the  other  m-sequence  at  non-peak  locations  (3.39),  its  extra  variance 
due  to  AWGN  from  N(0,  cr)  is  (cr6  +  3<J4  +  3<j2)/N  [1].  Although  not  yet  derived,  the 
corresponding  extra  variance  at  non-peak  locations  is  at  least  as  great.  Added  noise  clearly 
reduces  the  power  of  the  F-test  in  (3.51). 

As  an  example,  consider  the  case:  L  =  31,  A  =20,  SNR  =  0  (cr  =  1),  y=0.5.  The  separate 
variances  are 


<7 


2 

0 


N  v 


n-M 
L  , 


0.019 


[due  to  partial  sample], 


T 


t(»H 

L 


0.034  [extra due  to  non -peak], 


(3.52) 


cr6  +  3cr4  +  3cr2 
N 


0.35 


[AWGN], 


A  small  set  of  simulations  resulted  in  a  constituent  m-sequence  detection  rate  of  only  35%  at 
0  dB,  but  85%  at  3  dB,  though  these  rates  are  pessimistic.  20  TCF  values  were  used  in  both 
the  peak  set  and  control  set,  so  the  F-value  used  was  F19J9  =1.8.  A  reduction  in  the  95% 
confidence  level  of  the  test  would  improve  the  detection  rate  at  the  expense  of  more  false 
alarms.  A  larger  N  or  multiple  intercepts  would  improve  the  detection  rate. 

The  use  of  kurtosis  is  a  powerful  aid  in  detecting  mixtures  of  Gaussian  distributions.  Gold 
code  detection  depends  on  distinguishing  the  distributions  in  Figures  3.2  and  3.3:  the  locally 
most  powerful  test  that  a  sample  is  drawn  from  a  mixture  of  Gaussian  distributions  against 
the  hypothesis  that  it  is  drawn  from  a  single  Gaussian  population  is  based  on  kurtosis, 
jU4/ n\ ,  where  jUj  is  the  j- th  moment  of  the  sample  about  the  mean.  Such  a  test  is  being 
investigated.  Also,  as  the  non-peak  variance  (3.48,  3.49)  depends  on  n  and  L,  methods  to 
narrow  the  search  for  m-sequences  by  using  variance  information  to  limit  the  range  for  L  are 
also  under  investigation:  prior  knowledge  of  SNR  is  valuable. 

Once  a  constituent  m-sequence  has  been  detected,  a  search  is  carried  out  for  one  of  the  few 
other  m-sequences  that  may  complete  a  preferred  pair.  The  phases  of  the  m-sequences  do 
not  affect  the  above  detection  process.  However,  the  Gold  code  is  finally  identified  by 
maximizing  the  correlation  between  the  intercept  and  relatively  phase-shifted  versions  of  the 
detected  preferred  pair  of  m-sequences. 
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4. 


Higher-order  (>3)  Correlation  Functions  of  Gold  Codes 


The  set  of  L2  cyclic  rotations  of  all  L  of  the  Gold  codes  generated  by  a  preferred  pair  of 
m-sequences  of  length  L,  rPwk(0<j<L-\ ,  0<k<L-] ),  is  closed  with  respect  to  x{  1 ,  -1 }: 

=«  iVi+kuj+l]vi+n+c 

=  uiui+nVi+kVi+n+k 

=  “i+avi+j3  (4-1) 

=  wi~a 

1  +  J 

i.e.  wkw-+ll=w-+J 


It  follows  that,  in  general,  higher-order  correlation  functions  are  3-valued.  It  has  already  been 
shown  that  the  triple-correlation  function  is  3-valued  and  has  no  peaks  of  value  1  except  for  rare 
special  cases  (3.4.1).  For  order  4: 


Ck  {p, q,  r )  =  —  £  wf  „  wL,  wk 

L  i 


1 


+p  i+q  i+r 

P 


— — y  w 

L  j  l+nP  l+ncjr 


(4.2) 


For  order  5: 


CA'  (p,  q,  r,.v)  =  y£  wf  w-+p  w-+q  w-+l.  wf+s 

Lj  i 

=  —  S  wf+n  wi  wL 

+P  l+nqr  ,+* 

=  ~rlwfi+jwL 

Li  i 

1  V  m 

=  y  Svv,+A. 

L  i 


(4.3) 


While  the  above  correlations  are,  in  general,  3-valued,  a  4th  value  of  1  would  result  if,  for  example, 
P-  k  and  j  =  s  in  (4.3).  However,  such  values  are  impossible,  except  in  cases  when  L  =  3 N,  until 
9lh-order  correlation  functions  are  reached. 
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4.1  9th-order  Correlation  Functions  (NCFs) 


For  the  Gold  code  wk  (wf  =  utvj+k ),  let  (p,q)  and  (. r,s )  be  respective  TCF  peak  locations  for 
m-sequences  u  and  v.  Then  the  following  9th-order  correlation  function  (NCF)  values  are 
peaks: 


C k  {p,  q,  r,  p  +  r,  q  +  r,  s,  p  +  5,  q  +  5)  =  y  I  wk  wk+p  wk+q  wk+r  wk+p+r  wk+q+r  wk+s  wk+p+s  wk+q+x 

L,  i= 0 


(4.4) 


But  (p,q )  is  a  TCF  peak  location  for  u  and  not  for  v,  so 

wi  wi+pWi+q  ~  ^ i^i+k^ i+p^i+k+pM i+qV i+k+q 

=  {uiUi+PUi+q  )  Vi+kVi+k+pVi+k+q 
^i+k^i+k+p^i+k+q 

=  A+/>  j*  0 


Also 

Wj+r  ^;/+ p+r  W i+q+r  ~~  ^ i+r  ^i+k+r^  i+ p+r ^i+k+p+r1*  i+q+r  T i+k+q+r 

~  i^i+r^i+r+p^i+r+q  )  ^ i+k+r^ i+k+p+r^ 'i+k+q+r 
^/+A'+r  Vj+k+ p+r  ^  i+k+q+r 

=  vi+i+r  from  (4.5) 


and 


(4.5) 


(4.6) 


^’i+.v  ^’/+ p+s  W j+q+s 


~  M i+s^i+k+s^ i+p+s^i+k+p+s^i+q+s^i+k+q+s 
~  fai+s^ j+s+pM i+s+q  )  ^i+k+s ^i+k+p+s ^i+k+q+s 
~  ^i+k+s ^ i+k+ p+s^i+k+q+s 

=  vi+j+s  from  (4.5) 


(4.7) 


As  ( r,s )  is  a  TCF  peak  for  v,  (4.5),  (4.6)  and  (4.7)  result  in  the  following  expression  for  the 
NCF: 


,  ,  .  1 
C  {p,q,r,p  +  r,q  +  r,s,p  +  s,q  +  s)  =  —  Iv;xiv;j 

L  i= 0 

=  1 


-i-V;  . 


(4.8) 
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These  NCF  peaks  occur  for  any  (p,q )  TCF  peak  of  u  and  any  ( r,s )  TCF  peak  of  v  and  for  any 
permutation  of  p,  q,  r,p+r ,  q+r,  s,  p+s,  q+s  in  Ck( ).  All  permutations  yield  identical  values 
and  do  not  increase  the  power  of  tests  for  Gold  codes  constructed  from  u  and  v.  However, 
each  combination  of  the  (L-  l)/2  independent  u  peaks  and  the  (L-  l)/2  independent  v 
peaks  yield  independent  Gold  code  peaks,  a  total  of  (L  -  1  )V4.  Thus  if  (pi,qd, 
i  =  1,2,..,  (L-  l)/2,  and  (rj,Sj),j  =  1,2,..,  (L-  l)/2  are  the  respective  locations  of  u  and  v 
TCF-peaks,  each  of  the  (L  -  1  )2/4  NCF  locations: 

(pi ,  q< ,  rj ,  pi  +  rj ,  qi  +  i-j ,  Sj ,  pi  +  Sj ,  qi  +  sj),  1  <  ij  <  (L  -1  )/2  (4.9) 

is  a  peak.  Thus  the  NCF  is  4-valued.  The  peak-values  of  1  are  extremely  rare:  the  vast 
majority  of  the  L 9  NCF  locations  produce  one  of  the  three  values  d-\,  do  or  0\  (3.1),  as  in 
(4.1),  (4.2)  and  (4.3). 

Tests  for  Gold  codes  are  based  on  the  mean  of  the  (L  -  1  )2/4  values  at  the  NCF  locations 
above,  which  should  be  peak  values  of  1  if  a  Gold  code  is  present.  If  no  such  code  from  the 
( u,v )  family  is  present,  the  expected  value  of  the  mean  is  -1/L.  For  each  possible  sequence 
length  L,  a  search  may  be  made  for  each  preferred  pair  ( u,v ),  i.e.  each  set  G(u,v).  Each  Gold 
code  in  G(u,v)  has  the  same  NCF  peaks,  determined  by  (pi,qi)  and  (r„.v/),  so  the  search  space 
is  small. 


4.2  Statistics  of  the  NCF  of  Gold  codes  in  noise 

Gold  code  w  of  length  L  in  additive  Gaussian  noise  n  [~N(0,  a2)]  has  NCF  mean: 


^-9  (-/l »  ^2 » —  *  -/8  )  —  ^ 

i 

=  e| 


1  k 


I  (wi  +  ni  twi+n  +  n‘+h  )•  •  •  (w/+;8  +  ni+ 78 ) 


L,= i 

-£{  wtwl+J]  ...wi+J& 

L,i=l 


+n'W'+j  1  -Wi+j 8 
+ni+j1wiwi+j2...wi+j8 


+n  i+m  Wi  Wi+j\  •  •  •  Wi+J1 
+n  j  fij+  wi+  ...  Wj+ 


+n,nl+Jinl+J2wl+j3  ...w,+J8 


+nini+,  ...«/+»}] 


1  v 

7  Y,  wi H’ 
Li= i 


i+J  i 


(4.10) 
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as  all  other  terms  involve  at  least  one  nk  (and  no  for  77?  >  1)  for  keS={i,i+  j  i  +./g}  and 


Ek/(«/, nh  •  •  •  I  5  and  L  *k)\ 

=  E[nk]E[f(nhnl2..)]  =  0 

(4.11) 


Thus  the  following  (unbiased)  NCF-means  result  for  peak  and  non-peak  locations: 


|  1  if  0',, . 78)  isaNCFpeak 

j-l/L  if  ,  y8 )  is  not  a  NCF  peak 


If  (7.  J2, . .  .,/8)  is  a  NCF-peak  location: 


W-W-i  •  . 

rvl r  1+ J\  • 

~wi+js='i=>niwi+jr-w*n 

=  n,w, 

nl+nwiwl+h  ... 

■Wi+j 8  =  W/+/|Wi'+;i 

ni+jXWiWHj\  ••• 

M;'+77  —  ni+M  W'+78 

ninl+Jlwl+J2  ... 

w;+;8  =«,n,.+;-|w;w/+yi 

The  variance  of  the  NCF  of  w  +  77  at  these  peak  locations  is 

1  Ji 


var[C9(./'l,./'2,...,./8)]  =  E 


7 1 (w,-  +  nt  )(wi+jl  +  ni+j] ). . .  (w,+;g  +  77(+;8 )-  C9 
L  (=1 

2 


=  E 


1  \k 


U  l,=i 


I  (w,  +  "i  )k,+ h  +  ni+j  ). . .  (wi+j  +  ni+J  ) 


C  2 


=  -j2EliWi 


wi+h  . . . wi+h  +...  +  wLwL+h  . . . wL+m 


+  7*1*1  +---  +  nLWL 

+  nUjlw]+j]  +...  +  nL+hwL+Jl 


+  ,7l+78  Wl+78  +  •  •  •  +  nL+j%  WL+j% 

+  n\n\+n  Wl+Jl  +  •  •  •  +  nLnL+n  WLWL+j\ 
+  nlnl+J2  w{wl+j2  + . . .  +  nLnL+h  wLwL+h 


(4.13) 


+  n\nl -"n\ +j%  +  n2n2+  j]  ■"n2+j8  +•••  + 


nLnL+h  ••• nL+J 


-1 


(4.14) 
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Squaring  E[]  in  (4.14)  and  omitting  all  cross-product  terms  involving  single  n*  expressions 
for  any  k,  which  have  zero  expectation,  gives  the  following  variance: 


var 


[Q> 0]  >  Ji ’•  •  •  ’  h )] ~  j2  ^ t 


[L2] 


+  w,  w2  wl+h  w2+Jj  . . .  wl+m  w2+jh  +  w,  w3 .  [0] 


+  I n2 wf  +  Z n2+J]  w2+h  +  ...+  Z n~+j% w2+j% 

1=1  /=i  (=i 

[9Lcr2] 

+  in2n2+nw2w2+n 

1=1 

[La4] 

+  Z  n2n2+j2w2w2+h 
/= 1 

[La4] 

+ 

> 

+  Z  ni+j2  ni+j&  Wi+ji  wi+j% 

/'= 1 

[La4] 

+  &i?n2+j]n2+j2w2w2+jlw2+j2 

i= 1 

+ 


+  Z  ni+j(,  ni+jj  ni+j$  Wi+je  Wi+n  wi+j 8 

I~\ 

+ 


[La6] 


6!3! 


[La6] 


+  Zn?nlji 


/= i 


[Lcr18] 


-1 


91  ,91  * 

Lr  +9L(J"  +  -^-La4  +  —  La6  +...  +  La] 8 


1  I  t2  ,  nr  ,2 


2!7!  3!6! 


4s 


9! 


.2k 


L  k!(9-k)! 


(4.15) 


The  variance  of  the  NCF  at  peak  (and  also  non-peak)  locations  is  thus,  in  full: 
var  [C9(jt,j2,...,j8)\  = 

-{9a2  +  36a4  +  84a6  +  126a8  +  126a10  +  84a12  +  36a14  +  9a16+a18} 


(4.16) 
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5. 


Detection  Based  on  Triple  Products 


5.1  Detection  of  m-sequences 

An  m-sequence  u  may  be  detected  by  the  property  that,  if  (- p'-q )  is  a  TCF-peak  location, 
the  triple  product  'M,-  .  is  always  equal  to  one  (see  Figure  5.1): 

Vi.  ui=ui_p-ui_q>  =>  u jU j_p’Uj_(l'  —  1  (5.1) 


111  ...  1 

if  ip,q)  = 
{-p'-q) 


Q>',q), 

a  peak  for  u 


Figure  5.1  Triple-product  based  m-sequence  detector  (suggested  by  R  Gold) 

Each  m-sequence  (/  generator  polynomial)  of  length  L  may  be  represented  by  a  single  peak 
location,  and  each  location  examined  by  testing  the  hypothesis  that  the  outputs  in  Figure  5.1 
are  one. 

If  the  m-sequence  is  data  modulated  by  /?„  which  changes  every  Mth  value  (M  =  mV),  i.e. 

bi  =  ±  1,  jM+  1  <i<(j+  1)M,  0  <j  <  N  -  (5.2) 

then  the  output  of  the  detector  (Figure  5.1)  for  (p,q)  corresponding  to  a  TCF  w-peak,  and 
assuming  no  noise,  is  as  shown  in  Figure  5.2. 


DETECTOR 

delay  Qj,q) 

bibj-pbj-q 


Figure  5.2  Output  of  m-sequence  detector  after  data  modulation 
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^  3 

The  output  bjbj-pbj^  closely  resembles  the  data  signal  as  (+1)'  =  +1  and  (-1)‘  =  -1. 
Assuming  q>p,  errors  may  occur  in  relatively  few  b,  within  p  and  q  values  after  a  data 
boundary  where  the  previous  data  value  is  -1 ,  as  illustrated  in  Figure  5.3.  The 
representative  peak  should  be  chosen  to  ensure  q-p  is  as  small  as  possible:  although  q^p  (no 
such  peaks  exist),  it  is  possible  that  q-f. j  =  1.  E.g.,  for  L  =  31,  [45]  has  a  peak  (17,18)  and 
[75]  has  a  peak  (1 9,20),  for  both  of  which  q-p  =  1 .  Consequently,  the  notch  shown  in  the 
rectangular  pulse  in  Figure  5.3  would  be  only  one  chip  wide.  For  longer  m-sequences,  the 
effects  of  such  errors  would  become  negligible. 


Figure  5.3  Output  approximation  to  data 

The  detector  thus  also  acts  as  a  decoder  which  does  not  require  phase  acquisition.  Such  a 
self-reference  decoder  would  be  dependent  only  on  the  chosen  shift  pair  for  the  particular 
code  and  not  on  the  chip  rate.  It  would  also  possess  LPD/LPI  by  chip-rate  hopping  - 
pseudo-randomly  varying  the  chip  rate.  This  idea  has  not  been  developed  here  but  seems 
worthy  of  future  investigation. 
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5.1.1  Mean  and  Variance  of  m-Sequence  Detector  Output 

Assuming  addive  Gaussian  noise  n,  ~  N(0,  o2),  the  mean  detector  output  is 

E[  (bjUj  +ni)(bi_pui_p  +  ni_p\bi_qui_q  +  nj_q)  ] 

=  E[  bjbi_p bl_ljulul_pul_q 

+nibi_pui_pbj_qui_q  +  nj_pbjujbi_qui_q  + 

+  ninj_pbi_qui_q  +  nini_qbi_pui_p  +  ni_pni_qbiui 
+  nini_pni_q  ] 


b,b,-pb^q 


{{p,q) corresponds  to  a  peak} 


(5.3) 


Note  that,  in  the  case  of  no  data  modulation  (V7  •  bt  —  l),  the  mean  is  1 . 


As  cross  products  vanish,  the  variance  of  the  detector  output  is 


E  [nf+n? 


2  2  2  2  2  2 
+  77“  +  77  /7  +  77  /7  +  77 

i-p  *  ni-q  '  rii  ni~p  '  ,Li-q  '  ni-p 


2  2  2  2  n 

ni-q  +nini_pni_q] 


3o~  +  3o4  + 


(5.4) 


This  is  the  variance  of  each  output  value,  irrespective  of  whether  the  input 
m-sequence  is  modulated  or  corresponds  to  a  TCF  peak.  It  is  greater  than  the 
TCF  variance  in  [1 ,  pA.8]  by  a  factor  of  L,  as  expected,  as  the  TCF  is  the  mean  of  L 
triple  products. 
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5.2  Detection  of  Gold  Codes 


A  Gold  code  w  =  uv,  where  u  and  v  are  m-sequences  with  respective  TCF  peaks  (-p  '-q  ) 
and  (- r'-s  ),  may  be  detected  by  testing  if  the  outputs  of  the  2-stage  detector  in  Figure  5.2 
are  equal  to  one. 


Figure  5.4  2-stage  triple-product  based  Gold  code  detector  (suggested  by  R  Gold) 


u  and  v  may  be  found  by  a  2-dimensional  search  for  (p',q)  and  (r',s): 


(- p'-q )  a  peak  for  u  =>  Stage-1  output  UjVjUj-p 'V,-p 'U^{ 'V,^, '  =  (5.5) 


Group  closure  =>  v,  v/-p'V;-^'=  v,v;+*  =  v,+/,  as  (-p'-q)  not  a  peak  for  v.  (5.6) 


(-r'-s)  a  peak  for  v  =>  also  a  peak  for y-shifted  v,+/,  by  TCF  phase-invariance 

=>  Vi  •  Stage-2  output  v/+/V/+/_,-'V, •+/.*'=  1 .  (5.7) 
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If  the  Gold  code  is  data  modulated  with  as  in  (5.2),  then  the  noiseless  output  of  the  2- 
stage  detector  (Figure  5.4)  for  ( p,q )  and  (r,s)  values  corresponding  to  peaks  is  as  shown  in 
Figure  5.5. 


Figure  5.5  Output  of  Gold  code  detector  after  data  modulation 

The  output  bibj-pbi^bj-rbi^-rbi^^bj-sbj-p^bi-t,-*  closely  resembles  the  data  signal  b,  as 
(+1  )9  =  +1  and  (-1  )9  =  -1 .  Error  patterns  are  more  complicated  than  those  for  the 
m-sequence  detector  illustrated  in  Figure  5.3.  An  example  is  given  in  Figure  5.6  for 
preferred  pair  [45]  and  [75]  with  (p,q )  =  (17,18)  and  ( r,s )  =  (19,20).  Again,  the  effect  of 
these  errors  would  be  greatly  reduced  for  longer  sequences.  A  decoder  similar  to  that 
suggested  in  5.1  would  suffer  from  greatly  increased  noise. 


4  6  17  is  i9  20  31  36  38  48  49  50 51 


Figure  5.6  Output  of  2-stage  detector  as  approximation  to  data 


5.2.1  Mean  and  Variance  of  Gold  Code  2-stage  Detector  Output 

Assuming  noise  n,  ~  N(0,  a2),  the  mean  of  the  detector  output  <7,  is 


d,  =  E  few,-  +  n, -%-pWi-p  +«;-„)••• 

=  bjbj_pbi_qbi_rbi_p_l.bi_q_rbi_sbi_p_sbi_q_s 


For  no  data  modulation,  this  value  is  always  1 . 


The  variance  is  Lx  greater  than  that  of  the  NCF  variance  (4.1 6): 


var  [di  ]  =  9a2  +  36a4  +  84  cr6  +  126cr8  +  126cr10  +  84tr12  +  36cr14  +  9  <r16  +  <x18 

(5.9) 


5.2.2  Statistical  Tests  to  Detect  Gold  Codes 

Clearly,  the  high  variance  in  (5.9)  suggest  the  averaging  of  many  output  values  to 
test  possible  {p,q)  and  (r,s)  pairs.  If  data  modulation  is  present,  such  a  test  may  be 
based  on  output  spectra.  Assuming  a  chip  duration  %  data  modulation  interval  Mr 
and  bj  random,  the  expected  Fourier  transform  (FT)  is  Mr&incifMf),  ignoring  the 
effect  of  the  signal’s  error  notches.  The  power  spectrum  envelope  of  the  2-stage- 
detector  output  for  correct  (p,q )  and  ( r,s )  pairs  is  approximately  as  shown  in 
Figure  5.6. 


Figure  5.7  Output  power  spectrum  for  correct  (p,q),  (rjs) 

When  ( p,q )  and  (r,s)  are  not  peaks,  the  output  of  the  2-stage  detector  for 
unmodulated  input  is 


.k...k  ...k  ... k  ...k  ...it 

-P- 


k  ...k 


W  W  W-  W-  W-  W  W  W  W- 
r *  i  rvi-p  rv  / — lj  rvi—r  v i-p-r  rvi-q-r  i-s  ryi-p-s  v i-q-s  i-J 


W: 


(5.10) 


for  some  £  and  j  (0  <  £,j<L  -  1),  from  the  closure  law  (4.1),  i.e.  a  shifted  Gold 
code  from  the  set  G{u,v)  (3.7).  The  output  for  a  ^-modulated  Gold  code  is  thus 


(5.11) 


31 


which  may  be  approximated  by 


Mi 


wi-jbj  ■ 


(5.12) 


Now  the  FT  of  //'  is  given  by  the  convolution  of  FTs: 

FT [wf_jbi)=  Wftj  * M t s\nc(fM t),  (5.13) 

where  the  Gold  code  FT,  W/ij ,  has  zeros  at  klz,  k  =  ±1 ,  ±2,. . .  .  There  is  no  general 

expression  for  \Vf_j :  Gold  codes  of  the  same  length  have  different  FTs  and  power 

spectra.  However,  the  FT  in  (5.1 3)  has  zeros  at  A/ r  rather  than  the  /c/Mr for  peak 
( p,q ),  ( r,s )  values.  This  feature,  in  addition  to  the  different  expected  values  at  zero 
frequency,  is  the  basis  for  selecting  (p,q )  and  ( r,s )  corresponding  to  a  particular 
preferred  (m,v)  pair/G(«,v)  set.  Smoothing  the  detector’s  output,  to  reduce  the  effect 
of  error  notches,  would  reduce  noise  in  the  estimated  power  spectrum. 

For  correct  shifts,  the  output  may  be  approximated  by  a  random  signal  (±1), 
changing  at  Mrintervals,  with  power  spectrum  M2T  sine 2(fM t).  Incorrect  shifts 
produce  an  approximately  pseudo-random  (±1)  output,  changing  at  ^intervals,  with 
approximate  power  spectrum  T  sine This  is  illustrated  in  the  graphs  of  power 
spectra  following  listings  of  programs  Goldetmods  and  Goldetmodsn  in  the 
appendix.  Power  at  0  frequency  is  consistently  higher  for  correct  shifts 
corresponding  to  peaks,  even  at  -6  dB.  This  feature  offers  a  good  test  for  the 
presence  of  a  modulated  Gold  code.  Among  the  tests  under  investigation  is  one 
based  on  the  statistic 


2 >;//, 


/ 


(5.14) 


where  x,  are  power-spectrum  values  and  f,  are  corresponding  frequencies,  y has  high 
expected  values  for  correct  shifts. 
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6. 


Results 


6. 1  Tests  based  on  Gold  code  TCFs 

Tests  were  based  on  the  differences  in  the  statistical  moments  of  TCF  values  on  peak 
sets  (sets  of  TCF  peak  locations  for  either  u  or  v)  with  those  on  non-peak  sets,  for 
L  =  31.  Kurtosis  produced  very  poor  discrimination.  The  mean  alone  was  better: 
peak  sets  had  TCF  means  close  to  -1/31  =  -0.032,  whereas  non-peak  sets  had  means 
close  to  -0.075  due  to  larger  numbers  of -9/31  than  7/31.  Variance  was  always 
lower  for  peak  sets,  as  was  skewness,  again  due  to  the  imbalance  of -9/31  and  7/31 
values  in  non-peak  sets.  A  linear  combination  of  the  first  three  moments  produced 
the  detection/false-alarm  rates  in  Table  6.1,  using  single  31 -length  samples.  The 
detection  rates  may  be  improved  by  averaging  over  several  samples. 


Noise  level 

Detection  rate 

False-alarm  rate 

6dB 

98 

1 

3  dB 

67 

25 

OdB 

56 

34 

Table  6.1  TCF  moment-based  detection  and  false-alarm  %  for  Gold  codes 


6.2  Tests  based  on  Gold  code  NCFs 

Extensive  simulations  were  carried  out  to  determine  detection  and  false-alarm  rates 
for  31 -length  Gold  codes  with  varying  noise  levels.  The  preferred  pair  of 
m-sequences  used  were  w[45]  and  v[75],  as  defined  in  (3.10)  and  (3.12).  The  % 
detection  and  false  alarm  rates  in  Table  6.2  are  based  on  single  31-length  samples  of 
the  Gold  code  in  noise.  All  combinations  of  the  two  sets  of  m-sequence  TCF-peak 
locations  are  used  in  the  tests,  as  described  in  section  4.1 ,  producing  a  mean  of  225 
NCF  values.  This  mean  ju  is  compared  with  a  threshold  of  0.9,  chosen  because  of 
upward  bias  of  NCF  values  due  to  higher-order  statistical  properties  of  the  noise 
generator  used. 
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The  hypotheses  tested  were: 

Ho:  Gold  sequence  from  G(u,v )  not  present  (if  //  <  0.9) 

Hp  Gold  sequence  from  G(u,v )  present  (if  ju  >  0.9) 

The  %  detection/false-alarm  rates  are  derived  from  1000  simulations  at  each  noise 
level.  For  false-alarms  rates,  non-peak  locations  were  used,  the  recorded  rate  being 
the  %  of  mean  values  exceeding  0.9.  Note  that,  using  (4.16),  the  approximate 
variance  of  NCF  values  is  0.27,  consistent  with  observed  values  in  the  simulations. 


Noise  level 

Detection  rate 

False-alarm  rate 

3  dB 

100 

0 

0  dB 

99.3 

0.0 

-3  dB 

77.9 

26.2 

-6  dB 

58.3 

46.6 

Table  6.2  Gold  code  detection  and  false  alarm  %  based  on  NCFs  of  31-samples 


Further  simulations  were  carried  out  basing  detection  and  false  alarms  on  averages 
of  NCF  values  over  10  independent  31 -length  samples.  The  improved  %  rates  are 
shown  in  Table  6.3. 


Noise  level 

Detection  rate 

False-alarm  rate 

3  dB 

0 

0  dB 

0 

-3  dB 

100 

7 

-6  dB 

74 

43 

Table  6.3  Gold  code  detection  and  false-alarm  %  based  on  10-averages  of 
31 -sample  NCFs 

For  127-length  Gold  codes,  the  increased  length  and  extra  NCF  locations,  reducing 
variance  by  a  factor  of  72,  makes  detection  easier.  From  (4.16),  at  -3  dB  variance  is 
increased  by  a  factor  of  38.5,  so  100%  detection  and  0%  false-alarms  are  expected. 
At  -6  dB,  variance  is  increased  (4.1 6)  by  a  factor  of  3822,  which  means  53 
127-samples  must  be  averaged  to  achieve  the  0  dB  performance  for  31 -length 
sequences,  i.e.  99%  detection,  0%  false-alarms. 
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6.3  Tests  based  on  2-stage  triple  products  of  Gold  codes 

A  repeating  Gold  code  was  generated  from  the  same  preferred  pair  (, u,v ),  as  in  6.2 
with  independent  additive  noise.  A  similar  upward  bias  to  that  described  in  6.2  was 
noted  in  the  mean  of  the  2-stage  detector  output  of  5.2.  The  same  threshold  of  0.9 
was  therefore  selected.  Assuming  fjN  is  the  mean  output  over  N  successive  values, 
the  following  hypotheses  were  tested: 

Ho  :  ( p,q )  and  (r,s)  correct  (if  //v  >0.9) 

Hi  :  (p,q)  and  (r,s)  incorrect  (if  /JN  <0.9) 

The  detection  %  shown  in  Table  6.4  is  calculated  from  the  number  of  values  >0.9  in 
200  simulations  when  correct  shift  pairs  are  used;  the  false-alarm  %  is  from  the 
number  of  values  >0.9  when  the  shift  pairs  are  incorrect  for  the  input  Gold  code. 


N  =  number  of  outputs 

Detection 

False-alarm 

averaged 

rate 

rate 

500 

77.5 

35.5 

1000 

85.5 

30.0 

2000 

93.0 

22.5 

Table  6.4  Detection  and  false  alarm  %  from  output  of 
2-state  detector  for  0  dB  noise 


Clearly,  from  Table  6.4,  averages  of  many  output  values  are  necessary  to  achieve 
moderate  detection  rates,  even  at  0  dB.  As  output  variance  is  greater  than  that  of  the 
NCF  estimator  by  a  factor  of  L(L  -  1  )2/ 4,  for  L  =  3 1  approximately  7000  output 
values  must  be  averaged  to  match  the  NCF  detection  rate  of  99%  at  0  dB  given  in 
Table  6.2.  As  infeasibly  large  averages  would  be  required  at  higher  noise  levels,  no 
further  simulations  were  made. 
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7. 


Conclusions 


The  TCF  method  for  Gold  code  detection  in  3.4  has  the  advantage  of  low  complexity  and  is 
two-stage.  The  TCF  locations  for  possible  m-sequences  of  appropriate  length  are  examined 
separately.  When  it  is  decided  that  an  m-sequence  u  is  one  of  the  preferred  pair  from  which 
the  Gold  code  is  constructed,  the  search  for  a  preferred  complimentary  v  is  narrowed  down: 
it  is  not  necessary  to  simultaneously  find  u  and  v.  Phases  of  u  and  v  are  also  irrelevant, 
reducing  the  dimensionality  of  the  search.  However,  the  method  is  very  sensitive  to  noise 
and,  for  good  detection  rates,  many  L-samples  are  necessary  to  allow  TCF  averaging.  The 
noise  problem  reduces  for  longer  Gold  codes  as  each  TCF  value  results  from  larger  triple 
averages  and  more  TCF  values  are  available  from  each  L-sample. 

The  NCF  method  is  the  most  successful  overall.  The  lowest-order  correlation  function  with 
Gold  code  peaks  is  nine.  Although  very  noisy,  there  are  (L-l)2/4  NCF  peak  locations  to 
average  over,  and  each  NCF  value  is  an  average  of  L  ninth-order  products.  Thus  the 
computational  load  is  not  high  as  the  NCF  is  only  evaluated  at  a  small  number  of  locations 
for  each  possible  preferred-pair.  Even  for  L=31  at  0  dB,  99.3%  detection  of  Gold  codes  is 
possible,  with  0.0%  false- alarms,  from  a  single  3 l-sample,  and  increasing  L  rapidly 
improves  detection.  However,  detection  rates  fall  sharply  for  higher  noise  levels.  All 
preferred  pairs  ( u,v )  for  each  L  must  be  tried,  but  phases  of  u,v  and  the  Gold  code  are 
irrelevant. 

The  third  method,  a  2-stage  triple-correlator  detector,  has  the  virtue  of  simplicity  and  is 
easily  implementable  in  hardware.  Three  relatively  shifted  versions  of  the  input  are 
multiplied  to  produce  an  output  of  triple  products,  three  relatively  shifted  versions  of  which 
producing  the  output  of  the  2nd  stage.  As  single  triples  are  used,  the  output  is  extremely 
noisy  and  very  large  numbers  must  be  averaged  to  produce  moderate  detection  at  even  low 
noise  levels.  For  data-modulated  Gold  codes,  spectral  tests  for  codes  are  possible:  spectra  of 
outputs  for  correct  shift  pairs  are  easily  distinguishable  from  spectra  for  incorrect  shift  pairs 
from  their  increased  low-frequency  power.  In  principle,  the  2-stage  detector  may  also  be 
used  as  a  phase-independent  decoder,  but  the  noisy  output  would  restrict  its  use  to  long 
codes. 

The  NCF  and  correlator  detectors  in  sections  4  and  5  have  the  advantage  that  the  presence  of 
further  codes  from  the  set  G(u,v)  strengthen  statistical  tests,  whereas  the  statistical  properties 
of  the  TCF  discriminants  of  section  3  are  undermined. 
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Program  listing  of  Matlab  M-files. 


gold.m 

function  y  =  gold(N, start) 

%  This  function  generates  a  Gold  code  pseudo-random  binary  sequence 
%  of  length  N. 

%  Inputs,  length  of  sequence  start  and  mask 
%  Mask  is  a  two  row  vector  representing  the  tap  points  of  both 
%  generators  used  in  the  Gold  code 

%  start  is  a  vector  of  length  N  and  is  the  starting  state  of  one  of  the  generators 
%  Output  is  a  vector,  y  of  length  N  containing  the  Gold  code 

switch  N 
case  4 

mask  =  [0  0  1  1  ;1  0  0  1]; 
case  5 

mask  =  [0001  1;1  000  1]; 
otherwise 

error(This  length  prbs  not  supported') 
end 

a  =  ones(l,N); 
b  =  start; 
if(start  —0) 

error('Starting  state  must  not  be  all  zero') 
end 

for  k=l:2AN-l, 

zl  =  a  &  mask(l, :); 
z2  =  b  &  mask(2,  :); 
for  i  =  N:  -1 :  2, 
zl(i-l)  =  xor(zl(i),zl(i-l )); 
z2(i-l)  =  xor(z2(i),z2(i-l)); 
end 

a=  [zl  (1  ),a(  1  :(N- 1 ))];  %  shift  right  operator 
b  =  [z2(I),b(l  :(N- 1 ))]; 

yl  =  xor(a(l),b(l));  %  generate  Gold  code  here 

% 

if  b(l)==l 
bbb=-l ; 
else 
bbb=1; 
end 
% 

if  yl(l)==l,  %  convert  to  a  signal  voltage  (T  =  -IV,  'O'  =  +1 V) 

y00  =  -i; 

else 

y(k)  =  l ; 
end 
end 
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Goldet.m 


%  'Goldet' 

%  Gold's  2-stage  method  to  detect  Gold  codes 

% 

%  specify  delay  pairs  (p.q)  and  (r,s)-  possible  pairs  (17,18)  and  (19,20) 

% 

P-17; 

q=18; 

r=19; 

s=20; 

% 

%  Input  noise  factor  fn  and  number  of  3 1  -length  Gold  code  cycles  nc 

% 

fn=0; 

nc=500; 

% 

x=[l  11-1-11-1  -1  -1-11-1-11-1111  -1  -1  -1-1-1-11-11111  I]; 

% 

%  Generate  nc  cycles  of  31 -length  Gold  code  in  y 

% 

for  i=l  :nc 
ia=(i-l)*31; 
for  j=l  :31 

y(ia+j)=x(j)+randn*fn; 

end 

end 

%  Calculate  triple  product  for  shift  pair  (p,q) 

% 

ni=nc*31-30; 
for  i=l  :ni 

z(i)=y(i)*y(i+p)*y(i+q); 

end 

% 

%  Calculate  triple  product  of  1st  stage  output  for  shift  pair  (r,s) 

% 

ni=ni-30; 
for  i=l  :ni 

t(i)=z(i)*z(i+r)*z(i+s); 

end 

% 

for  i=  1 : 1 24 
c(i)=t(i); 
end 
disp(c) 

% 

%  Calculate  mean  and  variance  of  2nd  stage  output 

% 

s=0; 

for  i=l  :ni 
s=s+t(i); 
end 
s=s/ni; 

% 

ss=0; 
for  i=l:ni 
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d=t(i)-s; 

ss=ss+(d*d); 

end 

ss=ss/ni; 

% 

disp('No.  of  output  values  averaged') 
disp(ni) 

dispC'MEAN  and  VARIANCE') 

disp(s) 

disp(ss) 

% 


NCFGo  Id  peaks,  m 

%  NCFGoldpeaks.m 

%  9th-order-CF  of  Gold  codes  using  all  combinations  of  TCF  peak  locations:  detection  % 

% 

fn=l 

nii=1000 

nav=10 

% 

pq(U)=3; 

pq(l,2)=5; 

pq(2,l)=6; 

pq(2,2)=10; 

pq(3,I)=12; 

pq(3,2)=20; 

pq(4, 1  )=24; 

pq(4,2)=9; 

pq(5,l)=17; 

pq(5,2)=l  8; 

pq(6,l)=7; 

pq(6,2)=16; 

pq(7 ,1  )=1 4; 

pq(7,2)=l ; 

pq(8,l)=28; 

pq(8,2)=2; 

pq(9,l)=25; 

pq(9,2)=4; 

pq(  1 0,1  )=  1 9; 

pq(10,2)=8; 

pq(l  1,1)=1 1; 

pq(l  1,2)=23; 

pq(  1 2,1  )=22; 

pq(12,2)=15; 

pq(13,l)=13; 

pq(13,2)=30; 

pq(14,l)=26; 

pq(14,2)=29; 

pq(15,l)=21; 

pq(15,2)=27; 

% 

rs(l,l)=l; 
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rs(l,2)=12; 

rs(2,I)=2; 

rs(2,2)=24; 

rs(3,l)=4; 

rs(3,2)=17; 

rs(4,l)=8; 

rs(4,2)=3; 

rs(5,l)=16; 

rs(5,2)=6; 

rs(6,l  )=5; 

rs(6,2)=28; 

rs(7,l)=10; 

rs(7,2)=25; 

rs(8,l)=20; 

rs(8,2)=19; 

rs(9,l)=9; 

rs(9,2)=7; 

rs(10,l)=18; 

rs(10,2)=14; 

rs(  1 1 , 1  )=  1 1 ; 

rs(l  1,2)=30; 

rs(12,l)=22; 

rs(12,2)=29; 

rs(  13,1  )=]  3; 

rs(13,2)=27; 

rs(  ]  4, 1  )=26; 

rs(14,2)=23; 

rs(  15,1  )=2 1 ; 

rs(15,2)=15; 

% 

mns=0; 
vss=0; 
nod=0; 
for  iii=l  :nii 

ss=0; 

% 

for  ipq=  1:15 
for  irs=l:15 

%  Set  values  p  q  r  p+r  q+r  s  p+s  q+s 

il=pq(ipq,l); 

i2=pq(ipq,2); 

i3=rs(irs,l); 

i4=il+i3; 

i5=i2+i3; 

i6=rs(irs,2); 

i7=il+i6; 

i8=i2+i6; 

% 

x=[l  11-1-11-1  -1  -1-11-1-11-1111  -1  -1  -1-1-1-11-11111  1]; 

% 

%  set  x(32-62)  etc  to  same  values  -  sequence  periodic  -  after  adding  noise 
% 

for  i=l  :31 

x(i+31  )=x(i)+randn*fn; 
x(i+62)-x(i+31); 
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x(i+93)=x(i+31); 

end 

% 

%  calculate  9th-0-CF  c(i  1  ,i2,i3,i4,i5,i6,i7,i8) 

% 

% 

s=0; 

for  k=l:31 

y=x(k)*x(k+il)*x(k+i2)*x(k+i3)*x(k+i4)*x(k+i5)*x(k+i6)*x(k+i7)*x(k+i8); 

s=s+y; 

end 

si(ipq,irs)=s/31 ; 

ss=ss+s/31; 

end 

end 

% 

mn=ss/225; 
amn(iii)=mn; 
if  mn  >  0.9 
nod=nod+l; 
else 

nod=nod+0; 

end 

% 

mns=mns+mn; 

% 

%disp(si) 

% 

vs=0; 

foripq=l:15 

forirs=l:15 

vs=vs+(si(ipq,irs)-mn)A2; 

end 

end 

vs=vs/225; 

vss=vss+vs; 

end 

% 

mns=mns/10; 

vss=vss/100; 

%disp('MEAN  VAR’) 

%disp(mns) 

%disp(vss) 

% 

nod=100*nod/nii; 
disp('Detection  %') 
disp(nod) 

%disp(amn) 

dispC***************************************') 

anii=nii/nav; 

nod=0; 

% 

for  ii=l  :anii 
s=0; 

n=nav*(ii-l); 
for  i=n+l  :n+nav 
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s=s+amn(i); 

end 

s=s/nav; 
bmn(ii)=s; 
if  s  >  0.9 
nod=nod+l ; 
else 

nod=nod+0; 

end 

end 

% 

disp('Delection  %  based  on  nav-averages') 

nod=100*nod/anii; 

disp(nod) 

disp(bmn) 

% 


NCFGoldpeaksm.m 

%  NCFGoldpeaksm.m 

%  9th-order-CF  of  Gold  codes  using  all  combinations  of  TCF  peak  locations:  detection  % 
%  based  on  NCF  mean  for  all  expected  peak  locations. 

%  10  Gold  codes  used  from  same  set  G(u,v). 

% 

clear  all 
fn=0.5 
nii=100 
nav=10 

% 

N=5; 

L=2AN-1; 

startl  =[  1  1  1  1  l];%Start  state  of  Gold  generator 

start2=[  1  1110]; 

start3=[l  1  10  0]; 

start4=[l  1  0  0  0]; 

start5=[l  0  0  0  0]; 

start6=[  1  0  0  0  1]; 

start7=[  1  0  0  1  0]; 

start8=[1  00  1  1]; 

start9=[l  0  1  0  0]; 

start  10=[1  0  10  1]; 

% 

%  Generate  Gold  codes  1-10 

% 

gl=gold(N,startl); 
g2=gold(N,start2); 
g3=gold(N,start3); 
g4=gold(N,start4); 
g5=gold(N,start5); 
g6=gold(N,start6); 
g7=gold(N,start7); 
g8=gold(N,start8); 
g9=gold(N,start9); 
gl  0-gold(N, start!  0); 

% 

for  i=l  :L 
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Xx(i)=gl(i)+g2(i)+g3(i)+g4(i)+g5(i)+g6(i)+g7(i)+g8(i)+g9(i)+gl0(i); 

end 

% 

% 

%  Input  all  9-order  peak  locations  for  Gold  codes 

% 

pq(l,l)=3; 

pq(l,2)=5; 

pq(2,l)=6; 

pq(2,2)=10; 

pq(3,l)=12; 

pq(3,2)=20; 

pq(4,l  )=24; 

pq(4,2)=9; 

pq(5J)=l  7; 

pq(5,2)=l  8; 

pq(6,1  )=7; 

pq(6,2)=16; 

pq(7,l  )=1 4; 

pq(7,2)=l ; 

pq(8,l)=28; 

pq(8,2)=2; 

pq(9,l)=25; 

pq(9,2)=4; 

pq(  1 0, 1  )=1 9; 

pq(10,2)=8; 

pq(l  1,1)=  1 1; 

pq(ll,2)=23; 

pq(12,l)=22; 

pq(12,2)=15; 

pq(13,l)=13; 

pq(13,2)=30; 

pq(  1 4,1  )=26; 

pq(14,2)=29; 

pq(15,l)=21; 

Pq(  1 5,2)=27 ; 

% 

rs(l,1)=l; 

rs(l,2)=12; 

rs(2,l)=2; 

rs(2,2)=24; 

rs(3,l)=4; 

rs(3,2)=17; 

rs(4,l)=8; 

rs(4,2)=3; 

rs(5,l)=16; 

rs(5,2)=6; 

rs(6,l)=5; 

rs(6,2)=28; 

rs(7,l)=10; 

rs(7,2)=25; 

rs(8,l)=20; 

rs(8,2)—  1 9; 

rs(9,l)=9; 

rs(9,2)=7; 

rs(  1 0, 1  )=1 8; 
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rs(10,2)=14; 
rs(l  1 ,1  )=1 1 ; 
rs(l  1,2)=30; 
rs(  ]  2, 1  )=22; 
rs(12,2)=29; 
rs(]  3,1  )=1 3; 
rs(13,2)=27; 
rs(14,l)=26; 
rs(14,2)=23; 
rs(15,l)=21 ; 
rs(15,2)=15; 

% 

mns=0; 
vss=0; 
nod=0; 
for  iii=l  :nii 

ss=0; 

% 

for  ipq=  1:15 
for  irs=  1:15 

%  Set  values  p  q  r  p+r  q+r  s  p+s  q+s 

il=pq(ipq,l); 

i2=pq(ipq,2); 

i3=rs(irs,  1 ); 

i4=il+i3; 

i5=i2+i3; 

i6=rs(irs,2); 

i7=il+i6; 

i8=i2+i6; 

% 

%  set  x(32-62)  etc  to  same  values  -  sequence  periodic  -  after  adding  noise 

% 

for  i=l:31 

x(i+31)=xx(i)+randn*fn; 

x(i+62)=x(i+31); 

x(i+93)=x(i+31); 

end 


%  calculate  9th-0-CF  c(i  1  ,i2,i3,i4,i5,i6,i7,i8) 


% 

s=0; 

for  k=l:31 

y=x(k)*x(k+il)*x(k+i2)*x(k+i3)*x(k+i4)N:x(k+i5)*x(k+i6)*x(k+i7)*x(k+i8); 

s=s+y; 

end 

si(ipq,irs)=s/31 ; 
ss=ss+s/3 1 ; 
end 
end 

% 

mn=ss/225; 
amn(iii)=mn; 
if  mn  >  0.9 
nod=nod+l ; 
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else 

nod=nod+0; 

end 

% 

mns=mns+mn; 

% 

%disp(si) 

% 

vs=0; 

foripq=1:15 
for  irs=l :  1 5 

vs=vs+(si(ipq  ,irs)-mn)A2; 
end 
end 

vs=vs/225; 

vss=vss+vs; 

end 

% 

mns=mns/10; 

vss=vss/100; 

%disp('MEAN  VAR’) 

%disp(mns) 

%disp(vss) 

% 

nod=100*nod/nii; 
disp(’Detection  %') 
disp(nod) 

%disp(amn) 

disp( '***************************************') 

anii=nii/nav; 

nod=0; 

% 

for  ii=1  :anii 
s=0; 

n=nav*(ii-l); 
for  i=n+]  :n+nav 
s=s+amn(i); 
end 

s=s/nav; 
bmn(ii)=s; 
if  s  >  0.9 
nod=nod+l ; 
else 

nod=nod+0; 

end 

end 

% 

disp(’Detection  %  based  on  nav-averages') 

nod=l  00*nod/anii; 

disp(nod) 

disp(bmn) 

% 
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NCFGoldn  onp.  m 


%  NCFGoldnonp.m 

%  9th-order-CF  of  Gold  codes  using  combinations  of  TCF  non-peak  locations:  detection  % 

% 

fn=l 

nii=1000 

nav=10 

% 

pq(U)=3; 

pq(l,2)=5; 

pq(2,l)=6; 

pq(2,2)=10; 

pq(3,l)=12; 

pq(3,2)=20; 

pq(4,l)=24; 

pq(4,2)=9; 

pq(5,l)=17; 

pq(5,2)=l  8; 

pq(6,l)=7; 

pq(6,2)=16; 

Pq(7 , 1  )=  1 4; 

pq(7,2)=l ; 

pq(8,l)=28; 

pq(8,2)=2; 

pq(9,l)=25; 

pq(9,2)=4; 

pq(  1 0, 1  )=  19; 

pq(10,2)=8; 

pq(U,l)=ll; 

pq(l  1,2)=23; 

pq(12,l)=22; 

pq(l  2,2)=  15; 

pq(l  3,1  )=  1 3; 

pq(13,2)=30; 

pq(14,l)=26; 

pq(14,2)=29; 

pq(l  5,1  )=2 1 ; 

pq(15,2)=27; 

% 

rs(  1,1  )=1 ; 

rs(l,2)=12; 

rs(2,l)=2; 

rs(2,2)=24; 

rs(3,l)=4; 

rs(3,2)=17; 

rs(4,l)=8; 

rs(4,2)=3; 

rs(5, 1 )— 1 6; 

rs(5,2)=6; 

rs(6,l)=5; 

rs(6,2)=28; 

rs(7,l)=10; 

rs(7,2)=25; 

rs(8,l)=20; 

rs(8,2)=19; 
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rs(9,l)=9; 
rs(9,2)=7; 
rs(  1 0, 1  )=1 8; 
rs(10,2)=14; 
rs(  1 1 , 1 )=  1 1 ; 
rs(l  1,2)=30; 
rs(12,l)=22; 
rs(12,2)=29; 
rs(13,l)=13; 
rs(13,2)=27; 
rs(  1 4, 1  )=26; 
rs(14,2)=23; 
rs(15,l)=21; 
rs(  1 5,2)=1 5; 
% 

mns=0; 
vss=0; 
nod=0; 
for  iii=l  :nii 


ss=0; 

% 

for  ipq=l:15 
for  irs=l :  1 5 

%  Set  values  p  q  r  p+r  q+r  s  p+s  q+s 
il=pq(ipq,l); 
i2=pq(ipq,2)+l; 
if  i2>30 
i2=2; 
end 

i3=rs(irs,l); 
i4=il+i3; 
i5=i2+i3; 
i6=rs(irs,2)+ 1 ; 
if  i6>30 
i6=2; 
end 

i7=il+i6; 

i8=i2+i6; 

% 

x=[  1  11-1-11-1  -1  -1-11-1-11-1111  -1  -1  -1-1-1-11-11111  1]; 

% 

%  set  x(32-62)  etc  to  same  values  -  sequence  periodic  -  after  adding  noise 

% 

for  i=l:31 

x(i+31)=x(i)+randn*fn; 
x(i+62)=x(i+31); 
x(i+93)=x(i+3 1 ); 
end 

% 

%  calculate  9th-0-CF C(il,i2,i3,i4,i5,i6,i7,i8) 

% 

% 

s=0; 

for  k=  1:31 

y=x(k)*x(k+il)*x(k+i2)*x(k+i3)*x(k+i4)*x(k+i5)*x(k+i6)*x(k+i7)*x(k+i8); 
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s=s+y; 

end 

si(ipq.irs)=s/31 ; 
ss=ss+s/3 1 ; 
end 
end 

% 

mn=ss/225; 
amn(iii)=mn; 
if  mn  >  0.9 
nod=nod+l; 
else 

nod=nod+0; 

end 

% 

mns=mns+mn; 

% 

%disp(si) 

% 

vs=0; 

for  ipq=l :  1 5 
for  irs=l :  1 5 

vs=vs+(si(ipq,irs)-mn)A2; 

end 

end 

vs=vs/225; 

vss=vss+vs; 

end 

% 

mns=mns/10; 

vss=vss/100; 

%disp('MEAN  VAR’) 

%disp(mns) 

%disp(vss) 

% 

nod=100*nod/nii; 
disp('Detection  %') 
disp(nod) 

%disp(amn) 

dispC***************************************'’ 

anii=nii/nav; 

nod=0; 

% 

for  ii=l:anii 
s=0; 

n=nav*(ii-l ); 
for  i=n+l  :n+nav 
s=s+amn(i); 
end 

s=s/nav; 
bmn(ii)=s; 
if  s  >  0.9 
nod=nod+l ; 
else 

nod=nod+0; 

end 
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end 

% 

disp('Detection  %  based  on  nav-averages') 

nod=100*nod/anii; 

disp(nod) 

disp(bmn) 

% 


TCFGoldlO.m 

%  TCF  of  10  Gold  codes  from  the  same  u  &  v 

% 

clear  all 
xnfac=0; 

N=5; 

L=2AN-1; 

startl  =[1111  1];  %  Starting  state  of  Gold  generator 

start2  =[11110]; 

starts  =  [1  1  10  0]; 

starl4  =  [  1  10  00]; 

start5  =  [1  0  0  0  0]; 

start6  =  [10001]; 

starl7  =  [10010]; 

start8  =  [1  00  1  1]; 

start9  =  [  1  0  100]; 

startl  0  =  [1  0  10  1]; 

% 

%  Generate  Gold  codes  1-10 

% 

gl=gold(N, startl); 
g2=gold(N,start2); 
g3=gold(N,start3); 
g4=gold(N,start4); 
g5=gold(N,starl5); 
g6=gold(N,start6); 
g7=gold(N,start7); 
g8=gold(N,start8); 
g9=gold(N,start9); 
gl  0=gold(N, startl  0); 

% 

for  i=l  :L 

X(i)=gl(i)+g2(i)+g3(i)+g4(i)+g5(i)+g6(i)+g7(i)+g8(i)+g9(i)+gl0(i); 

end 

%  set  x(32-62)  to  same  values  -  sequence  periodic 
% 

for  i=l:31 
x(i+31)=x(i); 
end 

% 

%  add  noise 

% 

for  i=l:62 

x(i)=x(i)+xnfac*randn; 

end 

% 
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%  calculate  TCF  c(i,j) 

% 

for i=  1 :30 
for  j=l  :30 
s=0; 

for  k=l:31 

y=x(k)*x(k+i)*x(k+j); 

s=s+y; 

end 

c(i,j)=s; 

end 

end 

% 

%  Print  tcf  array,  10  cols  at  a  time 

% 

disp('Columns  1-10') 
disp(c(l  :30,1 :10» 
dispCColumns  11-20') 
disp(c(l:30,l  1:20)) 
dispCColumns  21-30') 
disp(c(  1:30,2 1:30)) 
figure(l); 
contourf(c); 
figure(2); 
surf(c); 


TCFGoldstatp.m 
%  TCFGoldstatp.m 

%  TCF  of  Gold  codes  with  mean,  var,  m(4)  and  kurtosis  for  PEAK  set 

% 

fn=  1 .0; 

disp('Noise  factor') 
disp(fn) 

% 

%  generate  Gold  code  from  preferred  pair  [45],  [75]. 

% 

x=[l  11-1-11-1  -1  -1-11-1-11-1111  -1  -1  -1-1-1-11-11111  1]; 

% 

%  add  Gaussian  noise  and  repeat  sequence. 

% 

disp('Number  of  iterations') 
nitt=l  00 

% 
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cm=0; 

ca2=0; 

ca4=0; 

ckurt=0; 

% 

pq(1,1)=3; 

pq(],2)=5; 

pq(2J)=6; 

pq(2,2)=10; 

pq(3,l)=12; 

pq(3,2)=20; 

pq(4,1)=24; 

pq(4,2)=9; 

pq(5,1)=17; 

pq(5,2)=18; 

pq(6,1)=7; 

pq(6,2)=16; 

pq(7,l)=14; 

pq(7,2)=l; 

pq(8,l)=28; 

pq(8,2)=2; 

pq(9,l)=25; 

pq(9,2)=4; 

pq(10.])=19; 

pq(l  0,2)=8; 

pq(l  1,1)=1 1; 
pq(l  1,2)=23; 

Pq(  12,1  )=22; 

pq(12,2)=15; 
pq(13,l)=13; 
pq(l  3,2)=30; 
pq(14J)=26; 
pq(14,2)=29; 
pq(15,1)=21; 
pq(15,2)=27; 
% 

rs(  1 , 1  )=  1 ; 

rs(  1 ,2)=1 2; 

rs(2,l)=2; 

rs(2,2)=24; 

rs(3,l  )=4; 

rs(3,2)=17; 

rs(4,l)=8; 

rs(4,2)=3; 

rs(5,l)=16; 

rs(5,2)=6; 

rs(6,l)=5; 

rs(6,2)=28; 

rs(7,l)=10; 

rs(7,2)=25; 

rs(8,l)=20; 

rs(8,2)=19; 

rs(9,l)=9; 

rs(9,2)=7; 

rs(10,l)=18; 

rs(l  0,2)— 14; 


rs(l  ],1)=1 1; 
rs(  1 1 ,2)=30; 
rs(12,l)=22; 
rs(12,2)=29; 
rs(l  3,1  )=1 3; 
rs(13,2)=27; 
rs(  1 4, 1  )=26; 
rs(14,2)=23; 
rs(15,l)=21; 
rs(15,2)=15; 

% 

for  nit=  1  :nitt 

% 

for  i=]  :31 

y(i)=x(i)+randn*fn; 

y(i+31)=y(i); 

end 

% 

%  calculate  TCF  c(i,j) 

% 

for  i=l  :30 
for  j= 1:30 
s=0; 

for  k=l:31 

z=y(k)*y(k+i)*y(k+j); 

s=s+z; 

end 

c(i,j)=s; 

end 

end 

% 

%  print  TCF  array,  10  cols  at  a  time 

% 

%disp('Columns  1-10') 

%disp(c(l:30,l:10» 

%disp('Columns  1 1-20') 

%disp(c(l:30, 11:20)) 

%disp('Columns  21-30') 

%disp(c(l:30, 21:30)) 

%disp( '****************************************') 
% 

%  Calculate  statistics 

% 

am=0; 
for  i=  1 : 1 5 
ill=pq(i,l); 
i  1 2=pq(i,2); 
i2I=rs(i,l); 
i22=rs(i,2); 

am=am+c(i  1 1  ,i  1 2)+c(i2 1  ,i22); 
end 

mean=am/30; 

%disp('Mean’) 

%disp(mean) 

cm=cm+mean; 

zm(nit)=mean; 
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% 

a2=0; 
a4=0; 
for  i=l :  15 
il  l=pq(i,l); 
il2=pq(i,2); 
i21=rs(i,l); 
i22=rs(i,2); 
b2=c(il  l,il2)-mean; 
b2=b2*b2; 
a2=a2+b2; 
a4=a4+b2*b2; 
b2=c(i21  ,i22)-mean; 
b2=b2*b2; 
a2=a2+b2; 
a4=a4+b2*b2; 
end 

a2=a2/30; 
a4=a4/30; 
kurt=0; 
if  a2>0. 1 
kurt=a4/(a2*a2); 
end 

%disp('Variance’) 

%disp(a2) 

%disp('4th  moment  about  mean') 

%disp(a4) 

%disp('Kurtosis') 

%disp(kurt) 

ca2=ca2+a2; 

ca4=ca4+a4; 

ckurt=ckurt+kurt; 

% 

z2(nit)=a2; 

zk(nit)=kurt; 

% 

end 

% 

cm=cm/nitt; 

ca2=ca2/nitt; 

ca4=ca4/nitt; 

ckurt=ckurl/nitt; 

disp('OVERALL  MEANS  of  mean,  var,  m(4),  kurtosis') 

disp(cm) 

disp(ca2) 

disp(ca4) 

disp(ckurt) 

% 

disp('OVERALL  VARS  OF  mean,  var,  kurtosis') 

vm=0; 

v2=0; 

vk=0; 

% 

for  i=l  mitt 
d=zm(i)-cm; 
vm=vm+d*d; 
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d=z2(i)-ca2; 

v2=v2+d*d; 

d=zk(i)-ckurt; 

vk=vk+d*d; 

end 

% 

vm=vm/nitt; 

v2=v2/nitt; 

vk=vk/nitt; 

disp(vm) 

disp(v2) 

disp(vk) 

% 


T CFG  old st  a  tnp.  m 

%  TCFGoldstatnp.m 

%  TCF  with  mean,  var,  m(4)  and  kurtosis  -  for  non-peak  set 

% 

fn=l  .0; 

disp('Noise  factor') 
disp(fn) 

% 

%  generate  Gold  code  from  preferred  pair  [45],  [75]. 

% 

x=[l  11-1-11-1  -1  -1-11-1-11-1111  -1  -1  -1-1-1-11-11111  1]; 

% 

%  add  Gaussian  noise  and  repeat  sequence. 

% 

disp('Number  of  iterations’) 
nitt=100 

% 

cm=0; 

ca2=0; 

ca4=0; 

ckurt=0; 

% 

pq(l,l)=3; 

pq(l,2)=5; 

pq(2,l)=6; 

Pq(2,2)=10; 

pq(3,l)=12; 

pq(3,2)=20; 

pq(4,l)=24; 

pq(4,2)=9; 

pq(5,l)=17; 

pq(5,2)=l  8; 

pq(6,l)=7; 

pq(6,2)=16; 

pq(7,l)=14; 

pq(7,2)=l; 

pq(8,l)=28; 

pq(8,2)=2; 

pq(9,l)=25; 

pq(9,2)=4; 
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pq(  1 0,1  )=1 9; 
pq(10,2)=8; 
pq(ll,l)=l  1; 
pq(l  1,2)=23; 
pq(  1 2, 1  )=22; 
pq(12,2)=15; 
pq(13,l)=13; 
pq(13,2)=30; 
pq(14,l)=26; 
pq(]4,2)=29; 
pq(]  5,1  )=21 ; 
pq(15,2)=27; 

% 

rs( 1,1 )= 1 ; 

rs(],2)=12; 

rs(2,l)=2; 

rs(2,2)=24; 

rs(3,l)=4; 

rs(3,2)=17; 

rs(4,l)=8; 

rs(4,2)=3; 

rs(5,l)=16; 

rs(5,2)=6; 

rs(6,l)=5; 

rs(6,2)=28; 

rs(7,l)=10; 

rs(7,2)=25; 

rs(8,l)=20; 

rs(8,2)=19; 

rs(9,l)=9; 

rs(9,2)=7; 

rs(  1 0, 1  )=1 8; 

rs(10,2)=14; 

rs(l  1,1)=]  1; 

rs(]  1 ,2)=30; 

rs(  1 2, 1  )=22; 

rs(12,2)=29; 

rs(  13,1  )=  1 3; 

rs(13,2)=27; 

rs(  ]  4, 1  )=26; 

rs(14,2)=23; 

rs(15,])=21 ; 

rs(15,2)=15; 

% 

fornit=l:nitt 

% 

for  i=l  :31 
y(i)=x(i)+randn*fn; 
y(i+3i)=y(i); 
end 

% 

%  calculate  TCF  c(ij) 

% 

for  i=l  :30 
for  j=l  :30 
s=0; 
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for  k=l  :31 

z=y(k)*y(k+i)*y(k+j); 

s=s+z; 

end 

c(ij)=s; 

end 

end 

% 

%  print  TCF  array,  10  cols  at  a  time 

% 

%disp('Columns  1-10’) 

%disp(c(l  :30, 1 : 1 0)) 

%disp('Columns  11-20') 

%disp(c(l:  30, 11:20)) 

%disp('Columns  21-30') 

%disp(c(l:30, 21:30)) 

%disp('****************************!lt  ***********') 

% 

%  Calculate  statistics 

% 

am=0; 
for  i=  1 : 1 5 
il  l=pq(i,l); 
il2=pq(i,2)+l; 
if  i  1 2>30; 

i  1 2=2; 
end 

i2i=rs(i,l); 
i22=rs(i,2)+l ; 
if  i22>30 
i22=2; 
end 

am=am+c(i  1 1  ,i  1 2)+c(i2 1  ,i22); 
end 

mean=am/30; 

%disp('Mean') 

%disp(mean) 

cm=cm+mean; 

zm(nit)=mean; 

% 

%  create  non-peak  location  -  replace  q  by  q+1  or  2  if  q+l>30;  same  for  s. 
a2=0; 
a4=0; 
for  i=  1 : 1 5 
ill=pq(i,l); 
il2=pq(i,2)+l; 
if  i  1 2>30 
i  12=2; 
end 

i2I=rs(i,l); 
i22=rs(i,2)+l ; 
if  i22>30 
i22=2; 
end 

b2=c(il  l,il2)-mcan; 
b2=b2*b2; 
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a2=a2+b2; 

a4=a4+b2*b2; 

b2=c(i21,i22)-mean; 

b2=b2*b2; 

a2=a2+b2; 

a4=a4+b2*b2; 

end 

a2=a2/30; 
a4=a4/30; 
kurt=0; 
if  a2>0. 1 
kurt=a4/(a2*a2); 
end 

%disp('Variancc') 

%disp(a2) 

%disp(l4th  moment  about  mean') 

%disp(a4) 

%disp('Kurtosis’) 

%disp(kurt) 

% 

z2(nit)=a2; 

zk(nil)=kurt; 

% 

ca2=ca2+a2; 

ca4=ca4+a4; 

ckurt=ckurt+kurt; 

% 

end 

% 

cm=cm/nitt; 

ca2=ca2/nitt; 

ca4=ca4/nitt; 

ckurt=ckurt/nitt; 

disp('OVERALL  MEANS  of  mean,  var,  m(4),  kurtosis') 

disp(cm) 

disp(ca2) 

disp(ca4) 

disp(ckurt) 

% 

dispfOVERALL  VARS  of  mean,  var,  kurtosis’) 

% 

vm=0; 

v2=0; 

vk=0; 

% 

for  i=l:nitt 
d=zm(i)-cm; 
vm=vm+d*d; 
d=z2(i)-ca2; 
v2=v2+d*d; 
d=zk(i)-ckurt; 
vk=vk+d*d; 
end 
% 

vm=vm/nitt; 

v2=v2/nitt; 
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vk=vk/nitl; 

% 

disp(vm) 

disp(v2) 

disp(vk) 

% 


Goldetmods.m 

%  'Goldetmods' 

%  Gold’s  2-stage  method  to  detect  Gold  codes,  with  spectral  test 
%  Program  produces  spectrum  plot  of  detector’s  output,  which  approximates 
%  the  data  signal  if  the  delay  pairs  are  correct  for  the  Gold  code. 

% 

%  specify  delay  pairs  (p,q)  and  (r,s) 

% 

P=17; 

q=18; 

r=19; 

s=20; 

% 

%  Input  noise  factor  fn  and  number  of  31 -length  Gold  code  cycles  nc 

% 

disp('Noise  factor') 

f n=  1.0 

nc=40; 

% 

x=[l  11-1-11-1  -1  -1-11-1-11-1111  -1  -1  -1-1-1-11-11111  1]; 

% 

%  Generate  nc  cycles  of  b-modulated  31 -length  Gold  code  in  y 

% 

for  i=l  :nc 
b=-l ; 
xr=rand; 
if  xr>  0.5 
b=l; 
end 

bb(i)=b; 
ia=(i-l)*31; 
for  j=l  :31 

y(ia+j)=b*x(j)+randn*fn; 

end 

end 

%  Calculate  triple  product  for  shift  pair  (p,q) 

% 

ni=nc*31-30; 
for  i=l  :ni 

z(i)=y(i)*y(i+p)*y(i+q); 

end 

% 

%  Calculate  triple  product  of  1st  stage  output  for  shift  pair  (r,s) 

% 

ni=ni-30; 
for  i=l  :ni 

t(i)=z(i):|:z(i+r)*z(i-t-s); 
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end 

% 

%for  i=l:124 
%  c(i)=t(i); 

%  bbb(i)=bh(i); 

%end 

%disp(c) 

%  d  i  s  p( *  *  * *  *  *  *  *  *  *  *  *  *  *  *  * *  **  * * **  *  **  *  *  *  *  * *  *  *  *  *  •) 

%disp(bbb) 

%d  i  Sp(,5i:  ^  ^  ^ 5b*  51=  ^  Hc?T:  5!=  =|c  =!=  M1  ^  5j=^c  *  *) 

% 

%  Calculate  mean  and  variance  of  2nd  stage  output 

% 

s=0; 

for  i=l:ni 
s=s+t(i); 
end 
s=s/ni; 

% 

ss=0; 
for  i=l  :ni 
d=t(i)-s; 
ss=ss+(d*d); 
end 

ss=ss/ni; 

% 

disp('No.  of  output  values  averaged') 
disp(ni) 

dispCMEAN  and  VARIANCE’) 

disp(s) 

disp(ss) 

% 

fy=fft(t,  1 024); 

Py=fy.*conj(fy)/l  024; 

% 

for  i=l:128 
j=(i-l)*4; 
cs=0; 
for  k=  1:4 
cs=cs+Py(j+k); 
end 

Ps(i)=cs/4; 

end 

% 

f=1000*(0: 127)/ 1024; 
plot(f,Ps(  1 : 1 28)) 

title('Spectrum  of  output  for  correct  shift-pairs.  0  dB  noise  level'); 

% 


Goldetmodsn.m 

%  'Goldetmodsn' 

%  Gold's  2-stage  method  to  detect  Gold  codes,  with  spectrum  plot 

% 

%  specify  delay  pairs  (p,q)  and  (r,s)  -  non-peaks 


23 


% 

p=17; 
q=19; 
r=20; 
s=2I ; 

% 

%  Input  noise  factor  fn  and  number  of  31 -length  Gold  code  cycles  nc 

% 

disp('Noise  factor') 

fn=1.0 

nc=40; 

% 

x=[l  11-1-11-1  -1  -1-11-1-11-1111  -1  -1  -1-1-1-11-11111  1]; 

% 

%  Generate  nc  cycles  of  b-modulated  31 -length  Gold  code  in  y 

% 

for  i=l  :nc 
b=-l; 
xr=rand; 
if  xr>  0.5 
b=l; 
end 

bb(i)=b; 
ia=(i-l)*31; 
for  j=l  :31 

y(ia+j)=b*x(j)+randn*fn; 

end 

end 

%  Calculate  triple  product  for  shift  pair  (p,q) 

% 

ni=nc*31-30; 
for  i=l:ni 

z(i)=y(i)*y(i+p)*y(i+q); 

end 

% 

%  Calculate  triple  product  of  1st  stage  output  for  shift  pair  (r,s) 

% 

ni=ni-30; 
for  i=l  :ni 

t(i)=z(i)*z(i+r)*z(i+s); 

end 

% 

%for  i=l:124 
%  c(i)=t(i); 

%  bbb(i)=bb(i); 

%end 

%disp(c) 

(T^dispC************************************1) 

%disp(bbb) 

%disp( '************************************') 

% 

%  Calculate  mean  and  variance  of  2nd  stage  output 

% 

s=0; 

for  i=l:ni 
s=s+t(i); 
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end 

s=s/ni; 

% 

ss=0; 
for  i=l:ni 
d=t(i)-s; 
ss=ss+(d*d); 
end 

ss=ss/ni; 

% 

disp('No.  of  output  values  averaged') 
disp(ni) 

dispf'MEAN  and  VARIANCE) 

disp(s) 

disp(ss) 

% 

fy=fft(t,  1 024); 

Py=fy.*conj(fy)/l  024; 

% 

for  i=l :  1 28 
j=(i-l)*4; 
cs=0; 
for  k=l  :4 
cs=cs+Py(j+k); 
end 

Ps(i)=cs/4; 

end 

% 

f=1 000:|:(0: 1 27)/l  024; 
plot(f,Ps(  1 : 1 28)) 

% 
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Graph  produced  by  Matlab  routine  Goldetmods.m. 


Spectrum  of  output  for  correct  shift-pairs.  6  dB  noise  level 
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Spectrum  of  output  for  incorrect  shift-pairs.  6  dB  noise  level 


